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ABSTRACT 

No mechanisms have hitherto been conclusively demonstrated to be responsible for initiating 
optically-luminous nuclear (Scyfcrt) activity in local disk galaxies. Only a small minority of such galax- 
ies are visibly disturbed in optical starlight, with the observed disturbances being at best marginally 
stronger than those found in matched samples of inactive galaxies. Here, we report the first system- 
atic study of an optically-selected sample of twenty-three active galaxies in atomic hydrogen (HI) gas, 
which is the most sensitive and enduring tracer known of tidal interactions. Eighteen of these galaxies 
are (generally) classified as Seyferts, with over half (and perhaps all) having [OIII] luminosities within 
two orders of magnitude of Quasi-Stellar Objects. Only ^28% of these Seyfert galaxies are visibly dis- 
turbed in optical DSS2 images. By contrast, ^94% of the same galaxies arc disturbed in HI, in nearly 
all cases not just spatially but also kinematically on galactic (> 20 kpc) scales. In at least ^67% and 
perhaps up to ~94% of cases, the observed HI disturbances can be traced to tidal interactions with 
neighboring galaxies detected also in HI. The majority of these neighboring galaxies have projected 
separations of < 100 kpc and differ in radial velocities by < 100 km s _1 from their respective Seyfert 
galaxies, and many have optical luminosities ranging from the Small to Large Magellanic Clouds. 
In a companion paper, we show that only ~15% of a matched control sample of inactive galaxies 
display comparable HI disturbances. Our results suggest that: i) most Seyfert galaxies (with high 
nuclear luminosities) have experienced tidal interactions in the recent past; ii) in most cases, these 
tidal interactions are responsible for initiating events that lead to their nuclear activity. 
Subject headings: galaxies: active — galaxies: Seyfert — galaxies: interactions — galaxies: structure 
— galaxies: ISM 



1. INTRODUCTION 

Active Galactic Nuclei (AGNs) are believed to be the 
luminous visible evidence for the vigorous accretion of 
gas onto a central SuperMassive Black Hole (SMBH). 
In the case of Quasi-Stellar Objects (QSOs), adopting 
the convenient (but arbitrary) convention that QSOs 
have absolute _B-band (nonstellar) nuclear magnitudes 
M B < -23 (Veron-Cetty & Veron, 1998), the required 
mass-accretion rate, assuming a 10% efficiency for the 
conversion of inf ailing matter into radiation (e.g., Nor- 
man & Silk, 1983), is at least ~0.1 M yr" 1 . Lower- 
luminosity AGNs have correspondingly smaller mass- 
accretion rates, with those described in this paper having 
mass-accretion rates ranging down to about two orders 
of magnitude lower (§2). If the AGNs studied here have 

Electronic address: ck2v@virginia.edu 
Electronic address: jlim@asiaa. sinica. cdu.tw 
Electronic address: ywtang@asiaa. sinica. edu.tw 
Electronic address: pho@asiaa.sinica.cdu.tw 



lifetimes of ~10 yrs, the total amount of matter required 
to feed their activity is therefore ~10 5 -10 7 M . The ori- 
gin or source of this material, and how this material is 
brought into the sphere of influence of the central SMBH, 
is one of the longstanding unresolved problems in AGN 
research. 

Many mechanisms for providing fuel to the central 
SMBHs of active galaxies have been proposed (e.g., for 
a recent brief review, see Martini, 2004). This is par- 
ticularly true at relatively low accretion rates, where 
stellar disruptions due to close encounters with the cen- 
tral SMBH, dynamical friction on molecular clouds, stel- 
lar winds from massive stars or mass loss from evolved 
stars (e.g., Norman & Scoville, 1988), or turbulence in 
the interstellar medium (e.g., review by Wada, 2004), 
are among the proposed hypotheses. At higher accretion 
rates, non-axisymmetric gravitational potentials such as 
galactic bars (e.g., review by Combes, 2004) or gravi- 
tational interactions with neighboring galaxies (e.g., re- 
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view by Schmitt, 2004) may be needed to drive suffi- 
cient quantities of gas inwards to fuel the AGN. At even 
higher accretion rates, corresponding perhaps to that in 
QSOs, cannibalisms of gas-rich galaxies (minor mergers) 
or mergers between two gas-rich galaxies (major mergers) 
may be necessary to carry sufficient fuel to the AGN. 

Because of their more accessible spatial scales, most 
observational studies to date have focussed on whether 
galactic bars or galaxy-galaxy interactions are responsi- 
ble for triggering AGNs (in disk galaxies). Nearly all 
such studies have involved optically-selected samples of 
Seyfert galaxies, which make up only a very small frac- 
tion (roughly 1%) of relatively bright field galaxies (e.g., 
Huchra & Burg, 1992; Stcpanian, 1994). Seyfert galaxies 
are preferentially found among early-tpe spiral galaxies, 
and exhibit very broad emission lines (from their nuclei) 
compared with their inactive counterparts. They are fur- 
ther divided into different classes depending on the rel- 
ative intensity of their broader (~1000-15000 km s _1 ) 
component, which is dominated by permitted lines, to 
their narrower component (~200-2000 km s _1 ), which 
is dominated by forbidden lines. At the two extremes lie 
Seyfert 1 galaxies, which have a relatively strong broad 
component, and Seyfert 2 galaxies, which have a de- 
tectable narrow component only (in unpolarizcd light). 

Searches for an excess of bars in Seyfert compared 
with inactive galaxies have yielded either null or at best 
weakly suggestive results. For example, neither Ho, 
Filippcnko, & Sargent (1997) nor Mulchaey & Regan 
(1997) found an excess of large-scale (> 1 kpc) bars be- 
tween Seyfert and inactive galaxies, but Knapcn, Shlos- 
man, & Peletier (2000) found a marginal excess of such 
bars amounting to 79% ± 8% in Seyfert compared with 
59% ± 9% in inactive galaxies. Both Regan & Mulchaey 
(1999) and Martini & Pogge (1999) found that only a 
small fraction of Seyfert galaxies in their sample exhibit 
small-scale (< 1 kpc) bars, and suggested that such bars 
cannot therefore be the primary mechanism for fueling 
Seyfert activity. On the other hand, Laine et al. (2002) 
found that Seyfert galaxies have an excess or at least 
equal fraction of bars at practically all length scales (from 
tens of parsecs to tens of kiloparsecs) compared with 
inactive galaxies; taken together, 73% ± 6% of Seyfert 
galaxies have bars compared with 50% ± 7% of inactive 
galaxies in their sample. Maia, Machado, & Willmer 
(2003) find bars among 28% of Seyfert compared with 
14% of inactive galaxies in their control sample. By con- 
trast, Hunt & Malkan (2004) find that Seyfert 1 galaxies 
exhibit practically no bars, whereas only a small fraction 
(15% ±5%) of Seyfert 2 and a similar fraction of inactive 
galaxies in their control sample exhibit bars. 

Rather than bars, both Regan & Mulchaey (1999) 
and Martini & Pogge (1999), as well as Pogge & Mar- 
tini (2002), found that the majority if not all of the 
Seyfert galaxies in their sample exhibit nuclear dust spi- 
rals. These dust spirals may trace shocks or turbulence 
in the central circumnuclear gas disk, and hence corre- 
spond to sites of angular momentum dissipation. Martini 
et al. (2003), however, found that nuclear dust spirals 
occur with comparable frequency in both active and in- 
active galaxies. They concluded that no universal fueling 
mechanism has yet been demonstrated to operate in disk 
galaxies at spatial scales larger than ^100 pc from the 
AGN. 



Searches for an excess of Seyfert galaxies among (ap- 
parently) interacting galaxies compared with matched 
samples of normal or apparently isolated spiral galax- 
ies have yielded similarly ambiguous if not contradictory 
results. Dahari (1985) find only a marginal excess of 
Seyfert galaxies among interacting spiral galaxies, but a 
clear excess when the sample is restricted to just strongly 
interacting spiral galaxies. Keel et al. (1985) also find 
Seyfert galaxies more frequently in interacting compared 
with isolated spiral galaxies, but by contrast with Dahari 
(1985) find a smaller excess in the most disrupted sys- 
tems. On the other hand, Bushouse (1986) find a smaller 
fraction of Seyfert galaxies among violently interacting 
than isolated spiral galaxies. 

Despite the conflicting results, examinations of opti- 
cal images do reveal that the fraction of Seyfert galaxies 
that are visibly (or apparently) interacting in the opti- 
cal is relatively small. Indeed, even in observations with 
the Hubble Space Telescope, the vast majority of local 
Seyfert galaxies (Malkan, Gorjian, & Tarn, 1998), as well 
as QSO host galaxies at redshifts as high as z « 0.3 (e.g., 
Bahcall et al., 1997; Boyce et al., 1998), do not appear 
to be visibly disturbed in the optical. Apart from their 
nuclear activity, local Seyfert galaxies appear to other- 
wise resemble the general population of inactive, prefer- 
entially early-type spiral, galaxies, de Robertis, Hayhoe 
& Yee (1998) find at best a marginal difference in light 
asymmetries between Seyfert and inactive galaxies; us- 
ing the same sample but a more sophisticated analysis, 
Virani, de Robertis, & VanDalfscn (2000) find no differ- 
ences. On the other hand, Hunt & Malkan (2004) find 
that 30% ± 9% of Seyfert 1 (and a comparable fraction 
of LINERs and HII) and 63% ± 7% of Seyfert 2 galaxies 
exhibit isophotal twists compared with just 6% ± 4% of 
a matched sample of inactive galaxies. 

With little direct evidence for (more frequent) galaxy- 
galaxy interactions, studies have turned to whether 
Seyfert galaxies more frequently exhibit projected neigh- 
boring galaxies by comparison with matched samples of 
inactive galaxies (e.g., summary in Schmitt, 2004). Such 
an excess has been found in some studies (e.g., Stauf- 
fcr, 1982; Dahari, 1984; Rafanelli et al., 1995; Dultzin- 
Hacyan et al., 1999; Koulouridis et al., 2006), but only 
marginally or not at all in others (e.g., Fucntcs- Williams 
& Stocke, 1988; de Robertis, Yee, & Hayhoe, 1998). Se- 
lecting Seyfert galaxies from the Markarian catalog (see 
§2 for a brief description of this catalog), MacKcnty 
(1989) find a larger fraction of Seyfert with projected 
neighboring galaxies compared with inactive galaxies ly- 
ing in the close vicinity, but not with randomly selected 
inactive galaxies selected also from the same catalog. 
Laurikaincn & Salo (1995) find that Seyfert galaxies ap- 
pear in paired systems no more frequently than inactive 
galaxies, but have on average more projected neighboring 
galaxies. Schmitt (2001) did not find any differences in 
the fraction of galaxies with genuine neighboring galaxies 
(i.e., with known redshifts) among Seyfert, LINERs, HII, 
transition, and absorption-line galaxies when considering 
only those with similar morphological types. 

In many of the post-1990's (but not all; c.f. Rafanelli 
et al., 1995) studies, the excess of Seyfert with projected 
neighboring galaxies (where found) is seen only around 
Seyfert 2 but not Seyfert 1 galaxies (Laurikainen & Salo, 
1995; Dultzin-Hacyan et al., 1999; Koulouridis et al., 
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2006), or alternatively that Seyfert 2 galaxies have on 
average a richer environment (dc Robcrtis, Yee, & Hay- 
hoe, 1998). This supports an earlier study by Petrosian 
(1982) that found a difference between the environments 
of Seyfert 1 and Seyfert 2 galaxies. Both Chatzichristou 
(2002) and Hunt & Malkan (2004) suggest that Seyfert 2 
galaxies arc generally at an earlier stage of interactions or 
dynamical disturbances compared with Seyfert 1 galax- 
ies. Furthermore, Storchi-Bergmann et al. (2001) find 
more mergers or close interactions among Seyfert 2 galax- 
ies with recent circumnuclear star formation than those 
without, and suggest an evolutionary sequence whereby 
the former are seen earlier after an interaction. 

Many of the abovementioned studies are now recog- 
nized to contain differences in their selection of Seyfert 
and/or comparison samples of inactive galaxies, as well 
as differences in their methodologies, that may have given 
rise to their disparate results. Although the more recent 
studies purport to carefully address these problems, they 
nevertheless reach entirely different conclusions. For ex- 
ample, selecting both Seyfert (Huchra & Burg, 1992) and 
matched inactive galaxies from the CfA Redshift Survey 
(Huchra et al., 1983), de Robertis, Yee, & Hayhoe (1998) 
find no difference in the fraction of galaxies with pro- 
jected neighboring galaxies in either samples, but do find 
(as mentioned earlier) that Seyfert 2 are more likely to 
have closer projected neighboring galaxies than Seyfert 1 
galaxies. Based on the sample of Seyfert, LINER, tran- 
sition, absorption-line, and HII galaxies selected from 
the Revised Shapley-Ames Catalog of Bright Galaxies 
(Sandage & Tammann, 1981) as classified spectroscop- 
ically by Ho, Filippenko, & Sargent (1997) (containing 
mostly low nuclear luminosity galaxies), Schmitt (2001) 
found the same percentage among all the different types 
of galaxies with projected neighboring galaxies. Select- 
ing from the catalog of Lipovetsky et al. (1988) contain- 
ing most if not all of the then known Scyfert galaxies 
(not, as stated by Dultzin-Hacyan et al., 1999, based on 
the Second Byurakan Survey), both Dultzin-Hacyan et 
al. (1999) and Koulouridis et al. (2006) find an excess 
of Seyfert 2 but not Seyfert 1 with projected neighbor- 
ing galaxies compared with matched samples of inactive 
galaxies. 

Failing to find a difference between the fraction of 
Seyfert and inactive galaxies with projected neighbor- 
ing galaxies, and with de Robertis, Hayhoe & Yee (1998) 
finding only a marginal difference in light asymmetries 
between Seyfert and inactive galaxies, de Robcrtis, Yee, 
& Hayhoe (1998) suggested that Seyfert galaxies may 
have cannibalised minor companion galaxies. To test this 
idea, Corbin (2000) looked for morphological differences 
in the light distribution between Seyferts and galaxies 
having other nuclear spectral types (absorption, H II, 
and LINER). Finding no such differences, nor any quali- 
tative evidence that minor mergers have occurred among 
any of the Seyfert galaxies or LINERs, Corbin (2000) 
concluded that either such events occurred > 1 Gyr ago 
and are essentially complete or not at all. 

In summary, most Seyfert galaxies do not appear to 
be visibly disturbed let alone interacting with neigh- 
boring galaxies. Even if more Seyferts than inactive 
galaxies have projected neighboring or genuine compan- 
ion galaxies, it is not clear what fraction are actually 
involved in tidal interactions. In this paper, we apply 



the lessons learnt from studies which demonstrate that 
neutral atomic hydrogen (HI) gas can reveal gravita- 
tional interactions between galaxies not visible in optical 
starlight. The best known example is the M81 group of 
galaxies, the closest galaxy group next to our own Local 
Group. Optical observations show no visible evidence for 
interactions between M81 and its nearest neighbours, the 
starburst galaxy M82 and the inactive galaxy NGC 3077. 
HI observations, on the other hand, reveal gaseous fila- 
ments connecting all three galaxies that trace tidal in- 
teractions between these galaxies (Yun, Ho, & Lo, 1994). 
The same situation has also been seen in a small sam- 
ple of relatively nearby galaxies hosting QSOs (Lim & 
Ho, 1999). Why are tidal interactions more easily seen 
in HI gas than optical starlight? In normal (i.e., undis- 
turbed) spiral galaxies, the HI gas disk usually extends 
nearly twice as far out as the stellar optical disk. The 
outskirts of the HI disk is not as strongly gravitation- 
ally bound and therefore more susceptible to external 
perturbations than the stellar optical disk; when per- 
turbed, the outer regions of the HI disk also take longer 
to dynamically relax. This makes HI a particularly sen- 
sitive and enduring tracer of gravitational interactions 
between galaxies. Comprising simultaneous imaging and 
spectroscopy, HI observations provide not just measure- 
ments of the gas kinematics but also accurate redshifts 
of surrounding gas-rich galaxies. These information, to- 
gether with the spatial structure of the tidal features, 
can be used to cleanly identify or more confidently pick 
out possible interacting neighboring galaxies. 

In this paper, we report the complete results of the 
first systematic imaging study of the HI gas in a sample 
of Seyfert galaxies, the results for a small subsample of 
which were reported by Greene, Lim, & Ho (2004). In 
an accompanying paper (Tang et al., 2007), we report 
the corresponding results for a matched control sample 
of inactive galaxies. Both the Seyfert and control sam- 
ple are restricted to a narrow range of distances, ensur- 
ing that both samples are observed with the same uni- 
form sensitivity and detection threshold. Prior to our 
study, HI imaging studies of Seyfert galaxies have usu- 
ally targeted individual objects of special interest. The 
available HI images were therefore highly heterogeneous 
in both spatial resolution and detection threshold (e.g., 
see review by Brinks & Mundell, 1996). In particular, 
not all the HI images have sufficient brightness sensitiv- 
ity to detect spatially extended emission characteristic 
of tidal features. Furthermore, it was not clear whether 
the observed galaxies are representative of the general 
population of Seyfert galaxies. 

This paper is organized as follows. In §2 we describe 
our selection of the Seyfert galaxy sample, and in §3 our 
observations and data reduction. In §4 we present the 
results for the individual galaxies, and compile statisti- 
cal trends for the ensemble sample. We encourage those 
readers interested in just the statistical results of our 
study to nevertheless glance at the HI images shown in 
Figures 1-36, which capture the essence of this work. 
In §5 we present our case that tidal interactions initi- 
ate events that lead to optically-luminous Seyfert activ- 
ity in a large fraction of local disk galaxies. In §6 we 
discuss the implications of our results for the nature of 
the interaction, gas infall timcscalcs to fuel the central 
SMBH, and lifetimes of optically-luminous Seyfert ac- 
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tivity. In §7 we provide a concise summary of our re- 
sults and conclusions. The catalog from which we draw 
our sample of active galaxies assumes a Hubble con- 
stant of H Q = 50 km s" 1 Mpc" 1 , and we retain this 
for ease of comparison with other active galaxies (in- 
cluding QSOs) in this catalog. In the rest of our pa- 
per, unless otherwise stated, we assume a Hubble con- 
stant H a = 67 km s _1 Mpc -1 (and £7 = 1) closer to the 
presently accepted value, and for consistency with the 
earlier results for a small subsample reported by Greene, 
Lim, & Ho (2004). 

2. SAMPLE SELECTION 

The study reported in this paper is part of a larger 
program to address the prevalence of galaxy-galaxy in- 
teractions among the more optically-luminous AGNs in 
the Local Universe. In a separate study, we are conduct- 
ing a HI imaging survey of all known galaxies hosting 
QSOs (i.e., with absolute nuclear B-band magnitudes 
M B < -23 for H a = 50 km s" 1 Mpc" 1 ) at declina- 
tions 8 > -40° and rcdshifts z < 0.07 (Lim & Ho, 
1999; Lim et al., 2001). Here, we report an extension 
of that study to lower luminosity Scyfcrt nuclei, osten- 
sibly (see clarification below) — 19 > Mb > —23 (for 
H Q = 50 km s- 1 Mpc" 1 ). 

Prior to the Sloan Digitized Sky Survey (the results 
from which had not yet been fully published circa 2000) , 
a large fraction (roughly half or more, depending on 
the redshift interval) of local Seyfert galaxies (exclud- 
ing QSOs) in the Northern Hemisphere were found in a 
few major optical surveys. The first such large-scale sur- 
vey was the Byurakan Spectral Sky Survey, the results 
from which are more popularly known as the Markar- 
ian catalog (Markarian et al., 1989). A good fraction 
of the Seyfert galaxies selected in many of the optical 
studies mentioned in §1 were drawn (at least in part) 
from or can be found in this catalog. A much smaller 
number of Seyfert galaxies found (or rediscovered) in the 
(first) CfA Redshift Survey comprise among the most 
optically-luminous galaxies at their respective redshifts 
in the Local Universe (Huchra et al., 1983). About half 
of these galaxies are contained in the Markarian catalog. 
A large number of relatively low-luminosity AGNs (re- 
ferred to as "dwarf Seyfert nuclei) in the nearest galaxies 
(spanning a much lower redshift range than the Markar- 
ian or CfA Seyfert catalogs) were selected spectroscop- 
ically by Ho, Filippcnko, & Sargent (1997). Each of 
these surveys contain their own peculiar selection effects, 
and furthermore adopt somewhat different criteria (e.g., 
line ratios) in classifying galaxies into different activity 
classes (e.g., Seyfert 1, Seyfert 2, LINER, and inactive 
galaxies). Together with those found in other surveys 
or small-scale studies, the currently known population of 
optically-selected Seyfert galaxies therefore comprise a 
heterogeneously selected sample. Nevertheless, these ob- 
jects comprise the general population of Seyfert galaxies 
from which inferences have been drawn about their over- 
all nature. 

After looking through the available surveys, it turned 
out that none individually suited our needs. For exam- 
ple, the number of objects per unit redshift interval in the 
Markarian or CfA Seyfert catalogs are not by themselves 
sufficiently large, whereas the number of objects per unit 
redshift interval in the Ho, Filippcnko, & Sargent (1997) 



catalog do not contain enough luminous Seyfert nuclei, 
for our purpose. The only catalog that fulfilled our needs 
(at the time) was that of Veron-Cetty & Veron (1998), 
which collects together all the known AGNs (selected 
in different ways by the different studies) at the time 
of compilation (updated approximately every two years 
since 1991). By its very nature, this catalog allows us to 
draw a sample of Seyfert galaxies that are as closely rep- 
resentative of the known (heterogeneously selected) pop- 
ulation as is possible based on the most widely accepted 
classifications for their activity class. If the same phe- 
nomenon (e.g., tidal interactions) can be widely found in 
such a heterogeneous sample, this would consitute pow- 
erful evidence for its broader significance among the gen- 
eral population. 

The Veron-Cetty & Veron (1998) catalog lists (among 
other relevant parameters) the absolute B-band magni- 
tude of each galaxy in a relatively small aperture centered 
on the nucleus. To compile a reasonably large number 
of objects with listed absolute B-band (nuclear) magni- 
tudes closely approaching QSOs, we had to reach out to 
a redshift of z « 0.015. To restrict the sample to es- 
sentially uniform distances, we selected all twenty-seven 
disk galaxies listed in the Veron-Cetty & Veron (1998) 
catalog at 0.015 < z < 0.017 and with -19 > M B > -23 
(for H Q — 50 km s _1 Mpc" 1 ) in the Northern Hemi- 
sphere (S > 0°). To fill in the alloted observing time 
in the year 2000, we added one more galaxy to the sam- 
ple, UGC 3157, which was selected from Veron-Cetty & 
Veron (2000) catalog using the same criteria as described 
above. With the selected redshift and absolute magni- 
tude range, the Veron-Cetty & Veron (1998) catalog con- 
tains only two definite elliptical galaxies, Perseus A and 
NGC 315, both of which are well known radio galax- 
ies. Imaging the HI emission (if any) of these galaxies is 
difficult in the presence of their strong radio continuum 
emission. The Veron-Cetty & Veron (1998) catalog also 
includes objects classified as LINERs (Low Ionization 
Nuclear Emission Line Regions, which show relatively 
strong lines of low ionization species such as [OI] and 
[SII], and which may comprise low-luminosity AGNs), 
and also those classified as H II (which show narrow 
Balmer and [OIII] lines characteristic of H II regions). 
We included these galaxies as they are sometimes clas- 
sified in reliable studies or other often used catalogs as 
AGNs; e.g., NGC 1167, which is classified as a LINER in 
Veron-Cetty & Veron (1998), is classified as a Seyfert 2 
in Ho, Filippcnko, & Sargent (1997). Later on, we will 
compile the statistical results of our study separately for 
the entire sample and just those classified as Seyferts in 
Veron-Cetty & Veron (1998) for the reader to compare. 

Table 1 contains nineteen galaxies with nuclear spec- 
tral types classified as Seyfert 1, Scyfcrt 2, or intermedi- 
ate Seyfert type (between S1.0 and S2.0), two with uncer- 
tain Seyfert classifications, five classified as LINERs, and 
two as H II. We retain their names as listed in the Veron- 
Cetty & Veron (1998) catalog, except for UGC 3995 
which is referred to in that catalog as UGC 3995A (in 
most literature sources, UGC 3995 A is used to refer 
to the neighboring galaxy of UGC 3995). Note that 
NGC 5548 is also known as Markarian 1509, NGC 7469 
as Markarian 1514, and NGC 7679 as Markarian 534. 
All the other listed parameters (coordinates, redshifts, 
and absolute B-band magnitudes) are also from Veron- 
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Cetty & Veron (1998) (those for UGC 3157 from Veron- 
Cetty & Veron, 2000), except for Hubble types and in- 
clinations which are taken from the Hyperlcda database 
(http://leda.univ-lyonl.fr; Paturel et al., 2003). Note 
that our selection is entirely blind and hence unlikely 
to be skewed in any particular way to the HI properties 
of the selected galaxies. 

The absolute B-band (nuclear) magnitudes listed in 
Veron-Cetty & Veron (1998) are taken from the liter- 
ature, and are measured in different ways by different 
observers (e.g., different aperture sizes). The listed val- 
ues may also include an unknown contribution from the 
host galaxy: this is usually not a problem in the case of 
QSOs, which (far) outshines their host galaxies, but be- 
comes increasingly problematic with decreasing nuclear 
luminosities. In addition, the nonstellar nuclear contin- 
uum is usually not directly viewed in Seyfcrt 2 nuclei 
(e.g., Mulchaey et al., 1994; Malkan, Gorjian, & Tarn, 
1998). Instead of the continuum, Mulchaey et al. (1994) 
have shown that the [OIII] emission of Seyfert nuclei 
is a reliable proxy for their optical luminosities. The 
[OIII] emission appears to be isotropic in Seyfert galax- 
ies and well correlated with the nuclear UV continuum in 
Scyfert 1 galaxies, and therefore a reliable measure too 
of the nuclear UV continuum in Seyfert 2 galaxies. 

We have searched the literature for measurements of 
the nuclear [OIII] emission from our sample galaxies. 
Such measurements are available for a total of fourteen 
of our sample galaxies in Dahari & de Robertis (1988) 
and Whittle (1992), comprising three Seyfert 1, three 
Seyfert 1.5, six Scyfert 1.9/2, one LINER,' and one H II 
galaxy; for the seven objects that Dahari & dc Rober- 
tis (1988) and Whittle (1992) have in common, their 
quoted [OIII] intensities typically agree to within a factor 
of ^2 (sec Table 1). The [OIII] luminosities of these four- 
teen galaxies arc quite evenly spread over the range from 
~10 40 2 -10 42 - 3 ergs s" 1 (for H = 50 km s" 1 Mpc" 1 ), 
even if the LINER and HII galaxies are excluded; their 
absolute (nuclear) B-band magnitudes as listed in Veron- 
Cetty & Veron (1998) span essentially the entire range of 
our sample, from -19.0 to -21.6. There is no obvious dif- 
ference between the range of [OIII] luminosities for the 
different Seyfert classes, with both the least and most 
[OIII] luminous objects classified as Seyfert 2. The [OIII] 
luminosities for three other LINERs in our sample can 
be found in Ho, Filippenko, & Sargent (1997), including 
NGC 1167 which is classified in that study as a Seyfert 2; 
these objects span a lower range in [OIII] luminosities 
of ~10 39:5 -10 40 ' 2 ergs s' 1 (for H D = 50 km s" 1 Mpc" 1 ). 
Indeed, the bulk of the galaxies (including Seyfcrts) in 
Ho, Filippenko, & Sargent (1997) have [OIII] luminosi- 
ties below ~10 40 ergs s _1 . Local QSOs (following the 
definition of Veron-Cetty & Veron, 1998) typically have 
[OIII] luminosities of > 10 42 ergs s _1 (e.g., see Kauff- 
mann et al., 2003, scaled to H a = 50 km s _1 Mpc^ 1 ). At 
least ~63% of the Scyfert galaxies in our sample therefore 
have optical nuclear luminosities spanning a range that 
is up to about two orders of magnitude below the lower 
limit for QSOs. The remainder span the same range of 
absolute B-band (nuclear) magnitude as listed in Veron- 
Cetty & Veron (1998), and likely so the same range in 
optical nuclear luminosities. 

Our sample contains essentially all of the reliably iden- 
tified Seyfert (disk) galaxies known in their rcdshift in- 



terval as of 1998. It includes all seven Seyfert galaxies in 
the CfA Seyfert catalog (Huchra et al., 1983) within our 
selected redshift interval, and thus comprises a complete 
subsample from this catalog. These galaxies are listed 
separately in Table 2. Fifteen of the twenty- four galax- 
ies in our sample can be found in the Markarian catalog 
(Markarian et al., 1989), including twelve of the twenty- 
one objects classified as Seyferts (under the current clas- 
sification for their activity classes). These galaxies also 
comprise a complete subsample from the Markarian cata- 
log within our selected rcdshift interval, and again listed 
separately in Table 2. As we shall see, our results are 
the same if restricted to just these two subsamples or for 
our entire sample, emphasizing the broader significance 
of our results to the general population. 

3. OBSERVATIONS AND DATA REDUCTION 

We observed our sample galaxies with the Very Large 
Array (VLA) 1 of the National Radio Astronomical Ob- 
servatory (NRAO). The observations were conducted 
over four sessions from 2000 to 2004, with a total alloted 
observing time of ^80 hours, and at night to avoid radio 
interference from the Sun. Relevant details of the ob- 
servations, including the flux, secondary (for amplitude 
and phase calibration), and bandpass calibrators used 
for each target object, are summarized in Table 3. The 
pointing (phase) center of the telescope was occasionally 
offset from the position of the target object to better 
surprcss strong continuum sources close to the edge of 
the field that can otherwise limit the dynamic range of 
our maps. We used the most compact configuration of 
the VLA (the D-array) that in the 21-cm line of HI gas 
provides a synthesized beam of ^60". The latter cor- 
responds to a spatial resolution of ^20 kpc at the dis- 
tance to our sample galaxies, which is well matched to 
the average optical sizes of these galaxies. Our study 
is therefore tailored towards the detection of extended 
emission on galactic scales characteristic of tidal features. 
All the galaxies in our sample were observed apart from 
NGC 266, a LINER, as it had previously been imaged in 
HI also in D-array by Simkin et al. (1987) (as reported 
by Hibbard et al., 2001). 

At 21 cm, the primary beam of the VLA is ^32' at full- 
width half-maximum, corresponding to a linear diameter 
of ^660 kpc. We configured the correlator to record sig- 
nals in orthogonal circular polarizations, and in sixty- 
four spectral channels over a bandwidth of 6.25 MHz 
and hence a channel separation of 97.7 kHz. This cor- 
responds to a velocity separation between channels of 
~21.0 km s _1 (the actual velocity width of each channel 
is ~25.2 km s _1 ), and a full range in velocity spanned 
by all channels of ^1350 km s _1 . For the assumed Cos- 
mology, the velocity range spanned by the correlator cor- 
responds to a redshift interval of Az w 0.005 (distance 
range of ~20 Mpc), compared with the redshift inter- 
val of Az 0.002 for our sample galaxies. The central 
channel is set to the systemic heliocentric velocity of the 
HI line if previously measured (usually with a single-dish 
telescope), or if not then at its reported optical redshift. 

We targeted a 5a detection threshold in HI column 
density of ^3 x 10 19 cm~ 2 per synthesized beam in a 
given channel with a velocity width ~25.2 km s _1 . Be- 
low roughly this column density, HI gas appears to be 
ionized by intergalactic UV radiation (Bland-Hawthorn, 
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Freeman, & Quinn, 1997). We emphasize that this 
does not necessarily imply that we detect all the HI gas 
present, as a given parcel of HI gas may have a larger 
velocity width than that spanned by just one channel (in 
such a case, the HI column density per channel is lower 
than that targeted). The HI column density thresh- 
old corresponds to a minimum detectable gas mass of 
^7 x 10 7 M Q per synthesized beam in one channel. The 
targeted detection threshold corresponds to a rms flux 
density of ^0.45 mJy beam -1 in one channel, requiring 
an integration time of ^2 hrs (observing time of ~2.5 hrs) 
per galaxy. 

We performed all the data reduction in the standard 
manner using the NRAO AIPS package. We inspected 
the data carefully for radio frequency interference (RFI) , 
which when present were excised from the data. Oc- 
casionally the RFI coincided with part of the HI line 
from the target galaxy; this did not pose a serious prob- 
lem when the RFI was intermittent, but when present 
throughout much if not all of the observation seriously 
affected the detectability of that portion of the HI line. 
As a consequence, we were forced to discard two of our 
sample galaxies, Markarian 352 and Markarian 1510 (one 
Scyfcrt and one H II galaxy) , both of which were badly 
affected by RFI. Following bandpass calibration, as well 
as time-dependent amplitude and phase calibration, we 
made continuum maps (i.e., combining the data from all 
the useable channels) for each target galaxy. To cor- 
rect for residual amplitude and phase errors, we self- 
calibrated the data for several iterations first in phase, 
followed by one or two iterations in both amplitude and 
phase when the continuum signal detectable over the 
primary beam was sufficiently strong. We then sub- 
tracted the continuum emission from the visibility data 
by interpolating between line-free channels (determined 
through trial and error) on either side of the HI line. 
Finally, we made maps of each channel, and corrected 
for the primary beam response of the antennas. We 
also combined the channel maps (without primary beam 
correction) to make maps in total intensity (zeroth mo- 
ment) and intensity-weighted mean velocity (first mo- 
ment), where we tried several different combinations of 
spectral smoothing and minimum signal cutoff level in 
each channel so as to suppress the noise and bring out 
faint diffuse features. 

4. RESULTS 

We omit from further consideration Markarian 352 
(Scyfcrt 1) and Markarian 1510 (H II), which as men- 
tioned in §3 are badly corrupted by RFI. Leaving 
aside these galaxies, only two galaxies, Markarian 359 
(Seyfert 1) and 2E 0507+1626 (Seyfert 1.5), were not de- 
tected in HI, as described in Greene, Lim, & Ho (2004). 
In the case of Markarian 359, a strong continuum source 
detected in the same field limited the accuracy of our 
passband calibration, and hence resulting map sensi- 
tivity. Unfortunately, this leaves only two remaining 
Seyfert Is in our sample. We also leave out NGC 266 
(LINER) , which although previously observed in D-array 
as mentioned in §3, has a shallower HI image than the 
galaxies described here. The HI image of this galaxy 
(Simkin et al., 1987; Hibbard et al., 2001) reveals that 
NGC 266 is interacting with the Seyfert 2 galaxy Markar- 
ian 348 (which has a listed redshift of z = 0.014 in Veron- 



Cetty & Veron, 1998, placing it just outside our selected 
redshift range). Including this result would only further 
strengthen the conclusions we draw from our study. 

The results for five galaxies (four Seyferts and one H II 
galaxy) have previously been reported by Greene, Lim, 
& Ho (2004), and are not repeated here (except for a 
short dicussion of Markarian 1157 for clarification). This 
leaves a remaining eighteen galaxies in our sample, the 
results for which are presented below. These twenty- 
three galaxies comprise our ensemble sample from which 
we draw our statistical results. The measured HI prop- 
erties of these galaxies (including all those reported in 
Greene, Lim, & Ho, 2004) are summarized in Table 4, 
where we list the integrated HI flux density S(HI) and 
associated uncertainty, the systemic HI velocity v sys (HI), 
the corresponding HI gas mass M(HI), whether they ap- 
pear disturbed in HI and hence the interacting group to 
which they belong based on an objective definition of the 
observed disturbances (see §4.1), and whether they are 
visibly disturbed based on images from the 2nd Digitized 
Sky Survey (DSS2) . The HI flux densities were computed 
from the total intensity (zeroth moment) maps, and con- 
verted to HI gas masses in the usual fashion assuming 
that the HI emission is optically thin. 

We reached our targeted detection threshold for all 
our sample galaxies except for NGC 1167 (LINER), 
UGC 3995 (Seyfert 2), and Markarian 341 (unclassified 
Seyfert type). These three galaxies were more badly af- 
fected by intermittent RFI than the remaining sample, 
and so more of their data had to be discarded. This re- 
sulted in a rms noise fluctuation that is a factor of <~1.1 
higher for Markarian 341, and a factor of ~1.3 higher 
for both NGC 1167 and UGC 3995; the degradation in 
sensitivity is not severe, although (perhaps just coinci- 
dentally) UGC 3995 and NGC 1167 are the only two 
galaxies in our sample for which we did not detect HI 
disturbances. We also detected in HI a number of galax- 
ies in the same field of view as our target galaxies. Those 
neighboring galaxies lying within the primary beam (i.e., 
within a radius of ^15'.8 about phase center) are listed in 
Table 5, and are referred to by their names from optical 
catalogs or, when we fail to find them listed in these cat- 
alogs, their positional coordinates with the prefix HI. We 
also list their projected separation from their correspond- 
ing active galaxies, optical redshifts where available, as 
well as their corrected apparent B-band magnitudes and 
Hubble types as cataloged in the Hyperleda database un- 
less otherwise noted. Their measured HI properties are 
summarized in Table 6 (including all those reported in 
Greene, Lim, & Ho, 2004). 

The HI moment and channel maps of all the galax- 
ies in our ensemble sample, apart from those previously 
reported by Greene, Lim, & Ho (2004), are shown in Fig- 
ures 1-36. The HI moment maps of neighboring galaxies 
listed in Table 5 not already shown in Figures 1-36 are 
shown in Figures 37-43. For display purposes (i.e., to 
preserve a uniform rms sensitivity throughout our maps) , 
the maps shown in the next section have not been cor- 
rected for the primary beam attenuation (which is in any 
case negligible over the field of view shown). The inte- 
grated HI flux densities and computed HI gas masses 
listed in Tables 4 and 6 have, of course, been corrected 
for the primary beam attenuation. We have overlaid each 
integrated HI intensity (zeroth moment) map on an op- 
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tical (either the red or blue filters) image from the DSS2. 
The latter provides a relatively uniform set of optical im- 
ages for all our sample galaxies, and was used to search 
for any disturbances in optical starlight. 

4.1. Individual Galaxies 

We present here the results for each galaxy, highlight- 
ing features of most relevance to this work. Based on 
these HI maps, we can separate the individual galaxies 
in our ensemble sample into the following three groups. 

Group I. Those that clearly show HI tidal features tracing 
interactions with neighboring galaxies. These fea- 
tures are in the form of tidal bridges that connect 
the two interacting galaxies, or an extension from 
one or both galaxies in the direction of the other 
(i.e., an incomplete tidal bridge). In addition, fea- 
tures in the form of tidal tails comprising a pro- 
trusion or curved extension from one galaxy on the 
side away from the other interacting galaxy are usu- 
ally seen. Thirteen of our sample galaxies fall into 
this group. 

Group II. Those that clearly show both spatial and kine- 
matic disturbances in HI, but which cannot be 
directly (but are probably linked) to interactions 
with neighboring galaxies where detectable. These 
disturbances are in general less prominent than the 
tidal features observed in the first group. Four of 
our sample galaxies fall into this group. 

Group III. Those that show marginal or no detectable HI dis- 
turbances. Six of our sample galaxies fall into this 
group, with at least two exhibiting possible distur- 
bances. 

4.1.1. Group I: HI tidal features tracing interactions 
with neighboring galaxies 

Twelve of our sample galaxies exhibit HI tidal features 
that clearly trace interactions with neighboring galaxies. 
In nearly all cases, counterparts of these HI tidal features 
also arc seen in the interacting neighboring galaxies. By 
contrast, only two of these galaxies exhibit optical distur- 
bances in DSS2 images that can be traced to interactions 
with neighboring galaxies. 

Four galaxies — Markarian 1, NGC 7679, NGC 7682, 
and ARK 539, all Seyfert 2s — exhibit a complete 
HI tidal bridge with a neighboring galaxy. NGC 7679 
and NGC 7682 are interacting with each other, and 
are the only interacting Seyfert pair in our sample. Of 
these, only NGC 7679 appears disturbed in the optical. 
Another eight galaxies — NGC 5548, Markarian 461, 
Markarian 1157 (reported in Greene, Lim, & Ho, 2004), 
NGC 841, Markarian 573, NGC 7591, NGC 7469 (re- 
ported in Greene, Lim, & Ho, 2004), Markarian 1040, 
arranged in descending order of the observed prominence 
of their tidal features — exhibit cither incomplete or rel- 
atively short tidal bridges, and usually also tidal tails. 
One is a Seyfert 1 (Markarian 1040), one a Seyfert 1.5 
(NGC 5548), four Seyfert 2s (Markarian 461, Markar- 
ian 1157, Markarian 573, and NGC 7469), one LINER 
(NGC 841), and one an unclassified Seyfert (NGC 7591). 
Of these, only one (Markarian 1040) exhibits optical 



tidal features. NGC 5548 does not appear to be opti- 
cally disturbed in DSS2 images, but in a deep ground- 
based optical image was found to exhibit a linear feature 
that we confirm to be a tidal feature in HI. NGC 7469 
does not exhibit any appreciable optical disturbances, 
but its neighboring galaxy IC 5238 exhibits prominent 
tidal features produced presumably by interactions with 
NGC 7469. 

To this group we add Markarian 341, which shows 
a tidal tail in the optical, and forms an optically- 
overlapping pair with its larger and even more optically- 
disturbed neighbour NGC 169. We cannot spatially sepa- 
rate these two galaxies in our HI observations, but detect 
extended tidal features from NGC 169 that we attribute 
to interactions with Markarian 341. Thus, altogether, 
Group I contains thirteen galaxies, of which only three 
exhibit detectable tidal features or disturbances in opti- 
cal starlight based on DSS2 images. 

Apart from extended tidal features, the HI disks of 
the majority of these galaxies have kinematics that 
deviate strongly from normal rotating disks. Detected 
here on spatial scales of ^20 kpc or larger, the observed 
kinematic perturbations correspond to disturbances on 
a global scale. 

Markarian 1 

Markarian 1 (Mrk 1) is a Sb galaxy hosting a Seyfert 2 
nucleus. A neighboring galaxy with comparable size and 
brightness, NGC 451, is visible about 1'.9 (40 kpc) to the 
south-east of Mrk 1, differing in optical systemic veloc- 
ity by only ~100 km s _1 (all optical systemic velocities 
from the NASA Extragalactic Database, NED). Neither 
galaxies, however, exhibit any visible disturbances in the 
optical, nor docs Mrk 1 exhibit any detectable bar even 
at the angular resolution of the HST (Malkan, Gorjian, 
& Tam, 1998). 

Mrk 1 has previously been imaged in HI with the Gi- 
ant Meterwave Radio Telescope (GMRT) by Omar et al. 
(2002) at a factor of ^2 higher angular resolution than 
that attained here. It is instructive to compare their re- 
sults (in particular, their Figs. 5 and 6) with ours (Fig. 1). 
Omar et al. (2002) find that the HI morphology and 
kinematics of Mrk 1 are disturbed, whereas those of the 
neighboring galaxy NGC 451 appear normal. No direct 
link, however, can be made between the observed distur- 
bances in Mrk 1 and tidal interactions with NGC 451. 

On the other hand, our HI image (Fig. 1) shows that 
Mrk 1 and NGC 451 shares a common gaseous envelope. 
An examination of the channel maps (Fig. 2) reveals a 
tidal bridge connecting the two galaxies, as well as tidal 
tails on opposing sides of the two galaxies. By contrast, 
the gaseous disturbances detected by Omar et al. (2002) 
cannot be clearly associated with tidal features. Our 
HI image is a factor of ^5.6 times more sensitive in flux 
density than that of Omar et al. (2002); furthermore, our 
synthesized beam is a factor of ~4 larger in area, result- 
ing in a factor of ^22 higher brightness sensitivity than 
the HI image of Omar et al. (2002) . The higher bright- 
ness sensitivity of our HI map allows us to detect more 
extended and dimmer emission characteristic of tidal fea- 
tures. 

Consistent with the higher angular-resolution but 
lower sensitivity image of Omar et al. (2002), our image 
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reveals that the spatial distribution and kinematics of 
the HI disk in Mrk 1 is strongly perturbed. Although 
the largest variation in velocity is along the optical 
major axis, the velocity does not change from rcdshift 
to blucshift along this axis as would be expected of 
a normal rotating disk. Like Omar ct al. (2002) we 
detect relatively weak HI absorption at velocities lower 
than that spanned by the HI emission (Fig. 2). These 
absorption features are most likely due to the HI clouds 
close to the central region of the galaxy which absorbed 
the continuum photons from the AGN. Note that these 
absorptions do not affect the moment maps, and are not 
responsible for the irregular kinematics observed. By 
contrast, the interacting neighboring galaxy NGC 451 
exhibits what appears to be a normal rotating disk, 
as confirmed in the higher angular-resolution image of 
Omar et al. (2002). This trend towards more disturbed 
HI kinematics in the active galaxy compared with its 
interacting neighboring inactive galaxy is seen in a 
number of interacting pairs in our sample. 

NGC 7679 and NGC 7682 

NGC 7679 and NGC 7682 are a closely-separated pair 
of Seyfert 2 galaxies with virtually identical optical sys- 
temic velocities. Pairs of Seyfert galaxies are compara- 
tively rare, and this is the only such example in our sam- 
ple. NGC 7679 is a barred SO-a galaxy and NGC 7682 a 
barred Sab galaxy, with NGC 7679 lying ~4f5 (-99 kpc) 
to the south-west of NGC 7682. NGC 7679 is clearly 
disturbed in the optical, exhibiting a particularly promi- 
nent extension to the east of the galaxy; i.e., in a di- 
rection closely towards NGC 7682. On the other hand, 
NGC 7682 appears to be normal in the optical, even at 
the angular resolution of the HST (Malkan, Gorjian, & 
Tarn, 1998). 

Like Mrk 1, our HI image of NGC 7679 and NGC 7682 
(Fig. 3) reveals that the two galaxies share a common 
gaseous envelope. The channel maps (Fig. 4) clearly 
show a tidal bridge connecting the two galaxies, as well 
as tidal tails on roughly opposing sides of these galaxies. 
The HI kinematics of both galaxies also are disturbed, 
with the isovelocity contours exhibiting a twist along 
their HI major axes. 

Ark 539 

ARK 539 is a spiral galaxy with an unclassified Hub- 
ble type (in HyperLcda) hosting a Seyfert 2 nucleus. A 
relatively faint galaxy, PGC 2370360, with no previously 
reported redshift lies — 3'A (74 kpc) to the south-east of 
ARK 539. There is no evidence for optical disturbances 
in either galaxy. 

Our HI image (Fig. 5) and channel maps (Fig. 6) reveal 
a curved tidal bridge connecting the southern or south- 
western side of ARK 539 with the south-western side of 
PGC 2370360. In addition, the HI gas on the eastern 
side of ARK 539 curls towards the south in the direction 
of PGC 2370360. The HI kinematic axis of ARK 539 
twists with radial distance from the center of the galaxy, 
following the same curl as the HI spatial distribution. 
PGC 2370360 has a narrow HI width suggestive of a face- 
on galaxy in accord with its optical image. Although 
therefore more difficult to judge whether the HI gas in 



this galaxy also is disturbed kincmatically, its HI gas 
does not appear to be as strongly disturbed spatially as 
ARK 539. 

We also detected HI gas in two other galaxies in the 
same field, PGC 2369294 and HI 1828+503 (Fig. 41). 
PGC 2369294 lies ~P5 (-33 kpc) to the south-east of 
HI 1828+503. The systemic velocities of PGC 2369294 
and HI 1828+503 differ by about -250 km s" 1 from 
those of ARK 539 and PGC 2370360, and these two pairs 
of galaxies do not appear to be interacting with each 
other. PGC 2369294 and HI 1828+503 share a common 
HI envelope, with at least PGC 2369294 exhibiting 
disturbed HI kinematics indicative of interactions with 
HI 1828+503. 

NGC 5548 

NGC 5548 is a SO-a ring galaxy hosting a Seyfert 1.5 
nucleus. In the DSS2 image, the optical morphology of 
this galaxy appears to be asymmetric, with an extended 
spiral arm emerging from the western side of the galaxy 
towards the south before winding north. Malkan, Gor- 
jian, & Tarn (1998) noted that NGC 5548 appears nor- 
mal in their HST image. In the deep ground-based op- 
tical image, however, Tyson et al. (1998) detected a lin- 
ear feature (length ~35 kpc) with a low surface bright- 
ness emerging from NGC 5548 towards the south-east. 
Tyson et al. (1998) suggest that this feature is a remnant 
tidal tail from a past interaction involving NGC 5548. 
NGC 5548 does not have a detectable bar even at the an- 
gular resolution of the HST (Pogge & Martini, 2002). A 
relatively dim inclined disk galaxy, PGC 1725892, which 
has a comparable optical systemic velocity as NGC 5548, 
lies 6!0 (-132 kpc) to the south-west of NGC 5548. This 
galaxy does not appear to be visibly disturbed in the 
optical. 

We detected in HI both NGC 5548 and PGC 1725892, 
as well as the faint galaxy SDSS J141824.74+250650.7 
lying -5.'9 (-130 kpc) to the south-east of NGC 5548 
(Figs. 7 and 8). The galaxy SDSS J141824.74+250650.7 
differs in systemic velocity by —200 km s _1 from 
NGC 5548, and does not appear to be optically dis- 
turbed. Of these three galaxies, the HI gas in 
NGC 5548 appears to be the most strongly dis- 
turbed. NGC 5548 exhibits a prominent HI exten- 
sion on the southern side of the galaxy that points 
east towards SDSS J141824.74+250650.7. We iden- 
tify this feature as an incomplete tidal bridge pro- 
duced by interactions with SDSS J141824.74+250650.7. 
This feature lies close to, but is spatially displaced 
from, the low-surface-brightness linear optical feature 
found by Tyson et al. (1998) as mentioned above. 
SDSS J141824.74+250650.7 exhibits a protuberance in 
HI on its north-western side, which we identify as the 
counterpart of the same incomplete tidal bridge from 
NGC 5548. NGC 5548 also exhibits a less prominent ex- 
tended HI feature on the north-western side of the galaxy 
that may be a tidal tail produced by interactions with 
SDSS J141824.74+250650.7. More likely, this feature is 
an incomplete tidal bridge produced by interactions with 
PGC 1725892, as PGC 1725892 exhibits a prominent HI 
extension on its north-eastern side that we identify as an 
incomplete tidal bridge produced by interactions with 
NGC 5548. 
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Tyson ct al. (1998) found that NGC 5548 has a de 
Vaucouleurs (i.e., elliptical galaxy like) radial surface 
brightness profile, and suggest that this galaxy may 
comprise a recent merger between two gas-rich galaxies. 
Our HI image, however, suggests that NGC 5548 is more 
likely a disk galaxy that has been strongly disturbed 
by recent tidal interactions with PGC 1725892 and 
SDSS J141824.74+250650.7. 

Markarian 461 

Markarian 461 (Mrk 461) is a Sab galaxy with a 
Seyfert 2 nucleus. No optical disturbances nor bars are 
apparent in this galaxy even at the angular resolution 
of the HST (Pogge & Martini, 2002). A relatively dim 
galaxy, PGC 2044413, differing in optical systemic ve- 
locity by -150 km s" 1 , lies just -l.'l (-22 kpc) to the 
south-west of Mrk 461. A brighter inclined disk galaxy, 
PGC 48870, can be seen at a larger angular separation of 
— 3f9 (—82 kpc) also to the south-west of Mrk 461; there 
are no reported redshift measurements for this galaxy. 
Neither PGC 2044413 nor PGC 48870 appear to be dis- 
turbed in the optical. 

We detected both Mrk 461 and PGC 48870 in HI 
(Fig. 9), but apparently not PGC 2044413 as can be seen 
mostly clearly in the channel maps (Fig. 10). (With a 
different systemic velocity, PGC 2044413 should appear 
in a different set of channels than Mrk 461.) The HI 
kinematic axis of Mrk 461 is oriented along its optical 
major axis. The HI major axis, however, is elongated 
approximately orthogonal to the HI kinematic axis, 
intersecting PGC 2044413 and also in a direction closely 
intersecting PGC 48870. The HI gas in PGC 48870 
also is disturbed, with the HI major and kinematic axis 
twisting away from the optical major axis on the western 
side of PGC 48870 towards Mrk 461. We therefore 
attribute the elongated HI feature from Mrk 461 towards 
PGC 48870 to an incomplete tidal bridge, and the elon- 
gated HI feature from Mrk 461 away from PGC 48870 to 
a tidal tail, produced by tidal interactions between these 
two galaxies. We cannot rule out, however, that these 
features also may have been produced by interactions 
between Mrk 461 and PGC 2044413, which lies closer in 
projection than PGC 48870. 

Markarian 1157 

A detailed description of Markarian 1157 (Mrk 1157), 
which does not appear to be disturbed in the optical, can 
be found in Greene, Lim, & Ho (2004). Not immediately 
obvious in the HI images of both Mrk 1157 and its 
neighboring galaxy PGC 5848 shown in Greene, Lim, 
& Ho (2004) is that the north-west to south-east HI 
extensions seen in both galaxies lie along a line joining 
these galaxies; i.e., PGC 5848 lies -10!4 (-202 kpc) to 
the south-east of Mrk 1157. These extensions are closely 
aligned with the major HI kinematic axes but not with 
the corresponding optical major axes of both galaxies, 
and trace an incomplete tidal bridge and opposing tidal 
tails. 

NGC 841 

NGC 841 is barred Sab galaxy hosting a LINER. A 



comparably bright galaxy, NGC 834, differing in opti- 
cal systemic velocity by only —50 km s _1 , lies — 10.'7 
(-209 kpc) to the north-west of NGC 841. Neither 
NGC 841 nor NGC 834 appear to be optically disturbed 
in the DSS2 image. 

Our HI image of NGC 841 (Fig. 11) reveals a rotating 
gas disk that appears to be normal on the south-eastern 
side, but which is disturbed both spatially and kine- 
matically on the north-western side; i.e., on the side 
facing NGC 841. The HI disk of NGC 834 is also 
disturbed, with the HI kinematic axis having a different 
orientation than the optical major axis of the galaxy. 
In addition, the HI disk of NGC 834 appears to be 
disturbed (elongated) on the side facing towards and 
away from NGC 841. All these HI disturbances can 
also be seen in the channel maps (Figs. 12 and 13). We 
attribute the disturbed HI features in both NGC 841 
and NGC 834 on their sides facing each other as an 
incomplete tidal bridge, and the disturbed HI feature in 
NGC 834 on the side away from NGC 841 as a tidal tail. 

Markarian 573 

Markarian 573 (Mrk 573) is a double-barred SO-a 
galaxy hosting a Scyfcrt 2 nucleus. The nuclear bar is 
clearly defined in the HST image of Pogge & Martini 
(2002). A relatively dim galaxy, PGC 1226059, with no 
previously reported redshift, lies — 3.'8 (—83 kpc) to the 
north-east of Mrk 573. Neither galaxies appear to be op- 
tically disturbed in the DSS2 image, nor Mrk 573 even 
at the angular resolution of the HST (Malkan, Gorjian, 
k, Tarn, 1998; Pogge & Martini, 2002). 

We detected both Mrk 573 and PGC 1226059 in HI 
(Figs. 14 and 15). The HI emission of PGC 1226059 is 
not spatially resolved, but reveals that this galaxy has 
a systemic velocity that differs by only —100 km s _1 
from that of Mrk 573. The HI gas in Mrk 573 is clearly 
disturbed spatially and kincmatically. The HI kinematic 
axis is better aligned with the optical major axis on the 
north-eastern side of center, but twists away from this 
axis on the south-western side. In addition, a prominent 
extension can be seen towards the south-eastern side of 
Mrk 573. These disturbances are located on the side 
of Mrk 573 opposite to PGC 1226059, and likely traces 
a tidal tail produced by interactions with its northern 
neighbour. 

NGC 7591 

NGC 7591 is a barred Sb galaxy with an uncertain 
Scyfcrt classification in Veron-Cetty & Veron (1998). 
This galaxy does not appear to be disturbed in the opti- 
cal. Kim ct al. (1995) classify NGC 7591 as a LINER. On 
the other hand, Veron, Goncalves, & Veron-Cetty (1997) 
show that NGC 7591 has a composite Seyfert and H II 
nuclear spectrum. As Seyfert galaxies quite commonly 
exhibit enchanced circumnuclcar star formation, we will 
assume here that NGC 7591 is indeed a Seyfert galaxy. 
An inclined disk galaxy, PGC 214933, with no previously 
reported redshift, lies — 1.'8 (—39 kpc) to the south-west 
of NGC 7591. Amram et al. (1994) suggested that this 
galaxy is warped based on their optical image. 

Thomas et al. (2004) have presented a relatively shal- 
low HI image of NGC 7591, taken at a factor of —2 high 
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angular resolution in C-array with the VLA. In their im- 
age, PGC 214933 was not detected in HI. By compari- 
son, we detected not just PGC 214933, but also a rel- 
atively dim galaxy lying ^3'.6 (^77 kpc) to the east of 
NGC 7591, hereafter referred to as HI 2318+065 (Figs. 16 
and 17). 

The HI disk of NGC 7591 is clearly asymmetric 
(Fig. 16). The main HI body is much more extended on 
the north-western side than on the south-eastern side. 
In addition, the HI gas exhibits a prominent extension 
on the south-eastern side reminiscent of a tidal feature. 
NGC 7591 and its south-western neighbour PGC 214933 
share a common HI envelope, which in the channel 
maps (Fig. 17) is more clearly revealed as a tidal bridge. 
NGC 7591 and its eastern neighbour HI 2318+065 also 
are linked in HI, which again in the channel maps is more 
clearly revealed as a tidal bridge. This is one of the few 
active galaxies in our sample, together with NGC 5548, 
which are interacting with multiple neighboring galaxies. 

NGC 7469 

A detailed description of NGC 7469 can be found in 
Greene, Lim, & Ho (2004). NGC 7469 is disturbed in HI 
by tidal interactions with IC 5283, which is disturbed 
both optically and in HI. 

Markarian 1040 

Markarian 1040 (Mrk 1040) is a highly inclined Sbc 
galaxy hosting a Seyfert 1 nucleus. It forms an optically- 
overlapping pair with PGC 212995 (UGC 1935), which 
lies just ~19" (-7 kpc) to the north of Mrk 1040. 
PGC 212995 has a comparable optical redshift to, but 
is both smaller and dimmer than, Mrk 1040. 

It is not clear whether Mrk 1040 is disturbed in the 
optical. Although not reportedly disturbed in the HST 
image of Malkan, Gorjian, & Tarn (1998), in the DSS2 
image Mrk 1040 exhibits relatively dim and highly ex- 
tended spiral-like arms on both the north-eastern and 
south-western sides. These features extend approxi- 
mately equidistance from center and appear to be sym- 
metric, and hence may just be ordinary spiral arms. On 
the other hand, these spiral-like arms may actually be 
tidal features dragged outwards in the plane of the galac- 
tic disk. Ward & Wilson (1978) find that PGC 212995 
also exhibits strong but narrower emission lines than 
Mrk 1040, and suggest that the activity in PGC 212995 
is caused by interactions with Mrk 1040. In addition, 
Amram et al. (1992) find that the optical velocity field 
of PGC 212995 is disturbed, and suggest that this dis- 
turbance is caused by interactions with Mrk 1040. By 
contrast, Amram et al. (1992) attribute the asymmetric 
optical velocity field in Mrk 1040 to a large-scale and 
massive bar. 

Mrk 1040 and PGC 212995 cannot be spatially sep- 
arated at the angular resolution of our HI observation. 
In our HI map (Fig. 18), the centroid in integrated HI 
intensity as well as body of HI as a whole is centered 
on Mrk 1040, the larger galaxy. There is no evidence 
from our channel maps (Fig. 19) for any detectable 
emission from PGC 212995. The HI disk of Mrk 1040 
is clearly disturbed both spatially and kinematically, 
although these disturbances can be much more clearly 



seen in the channel rather than moment maps. A 
prominent HI feature can be seen on the eastern side 
of the galaxy, located beyond the extended sprial-like 
optical arm on this side of the galaxy. This HI feature 
curls north at its easternmost periphery; i.e., to the side 
where PGC 212995 is located. A similarly extended HI 
feature can be seen on the western side of the galaxy, 
again located beyond the extended sprial-like optical 
arm on this side of the galaxy. We attribute these 
features to tidal disturbances caused by interactions 
with PGC 212995. 

Markarian 341 

Markarian 341 (Mrk 341) is a SO-a galaxy that forms 
an optically-overlapping pair with NGC 169, a bigger 
and brighter galaxy located just ^23" (~7 kpc) to the 
north. These two galaxies have nearly identical optical 
rcdshiftes. NGC 169 is a Sab galaxy, and is clearly dis- 
turbed in the optical. It exhibits a prominent extension 
to the east that was presumably produced by tidal in- 
teractions with Mrk 341. Mrk 341 itself also exhibits 
a relatively dim extension to the south, which is on the 
side opposite to NGC 169 and hence a tidal tail. Mrk 341 
has an unclassified nuclear spectral type in Veron-Cetty 
& Veron (1998). Keel et al. (1985) suggest that Mrk 341 
hosts a LINER, whereas Veron, Goncalves, & Veron- 
Cetty (1997) find a composite H II and Seyfert nuclear 
spectra type, as in the case of NGC 7591. We shall hence- 
forth assume that Mrk 341, just like NGC 7591, is indeed 
a Seyfert galaxy. 

Mrk 341 and NGC 169 cannot be spatially separated 
at the angular resolution of our HI observation. In our 
HI map (Fig. 20), the centroid in integrated HI intensity 
as well as body of HI as a whole is centered on NGC 169, 
the larger galaxy. There is no evidence from our channel 
maps (Fig. 21) for any detectable emission from Mrk 341. 
The HI disk of NGC 169 is clearly disturbed both spa- 
tially and kinematically It shows a prominent extension 
to the east coinciding with the optical extension men- 
tioned above but extending to much larger distances. In 
addition, the HI kinematic axis twists away from the HI 
and optical major axes on both the eastern and western 
sides of the galaxy, and on both sides curl in a direction 
towards Mrk 341. All these disturbances can be most 
obviously linked to tidal interactions with Mrk 341. 

4.1.2. Group II: HI disturbances likely produced by 
interactions 

Four of our sample galaxies clearly exhibit both spa- 
tial and kinematic perturbations in HI, but unlike the 
situation in Group I their observed perturbations cannot 
be immediately linked to interactions with neighboring 
galaxies (if any). By contrast with their appearcnecs 
in HI, only one appears to be disturbed optically in 
DSS2 images. Two — MS 04595+0327 (Seyfert 1) and 
NGC 513 (Seyfert 1.9), both reported in Greene, Lim, 
& Ho (2004) — have neighboring galaxies with projec- 
tion separations of only about 50 kpc detected also in 
HI, providing circumstantial evidence that their distur- 
bances are indeed caused by tidal interactions with these 
galaxies. Neither, however, appear disturbed in the op- 
tical. One — UGC 1395 (Seyfert 1.9) — is connected 
to a separate HI feature that may be related to a pro- 
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jcctcd neighboring galaxy. The optical disk of this galaxy 
is asymmetric, suggesting a disturbance. The last - 
Markarian 1158 (H II galaxy), reported in Greene, Lim, 
& Ho (2004) — does not have any neighboring galax- 
ies detected in HI, and does not appear to be disturbed 
in the optical although clearly disturbed in HI. Although 
the observed perturbations in all four galaxies are in gen- 
eral less prominent than the tidal features observed in 
Group I, these perturbations nevertheless correspond to 
disturbances on a global scale. 

Detailed descriptions of MS 04595+0327, NGC 513, 
and Markarian 1158 can be found in Greene, Lim, & Ho 
(2004). Here, we present only the results for UGC 1395. 

UGC 1395 

UGC 1395 is a Sb ring galaxy hosting a Seyfert 1.9 
nucleus. Although not classified as a barred galaxy in 
Hyperlcda, McLeod & Rieke (1995) find a strong bar 
in the near-IR. UGC 1395 exhibits a relatively dim but 
highly extended spiral arm on its western side that gives 
this galaxy an asymmetric appearence in the optical. A 
relatively faint galaxy, PGC 1304037, with an unknown 
optical redshift lies -2.'9 (-65 kpc) west of UGC 1395; 
i.e., on the same side as the extended western spiral arm. 

Our HI image of UGC 1395 (Fig. 22) reveals that 
this galaxy is disturbed both spatially and kincmatically. 
First, the HI kinematic axis exhibits a prominent twist 
on the northern side of the galaxy. Second, UGC 1395 
is linked to a HI condensation lying — 2.'1 (—48 kpc) to 
the west. This HI condensation spans a different velocity 
range than UGC 1395; furthermore, an examination of 
the channel maps (Fig. 23) reveals that this condensa- 
tion is indeed a separate object. This HI condensation 
lies in close proximity to, but its centroid does not co- 
incide with, PGC 1304037; we believe this HI condensa- 
tion is most likely associated with PGC 1304037, but to 
be conservative we consider the possibility that it is not 
and refer to this condensation as HI 0155+066. Both the 
moment and channel (Fig. 23) maps reveal a tidal bridge 
between UGC 1395 and HI 0155+066. 

4.1.3. Group III: Weak or no HI disturbances 

Three of the galaxies in our sample — Markarian 993 
(Seyfert 1.5), UGC 3157 (Seyfert 2) and Markarian 1419 
(LINER) — exhibit relatively weak HI disturbances. 
Two of these (UGC 3157 and Markarian 1419) appear 
to be disturbed optically in DSS2 images, with the 
possible HI disturbances apparently coinciding with 
the optical disturbances. Two (Markarian 993 and 
UGC 3157) have neighboring galaxies detected in HI 
that are excellent candidates for having produced the 
observed disturbances in the active galaxies. Another 
galaxy in our sample — IRAS 14082+1347 (LINER) 
- also appears to exhibit HI disturbances, but this 
galaxy was detected at a relatively low signal-to-noise 
ratio. It does not appear to be optically disturbed in 
DSS2 images. The remaining two galaxies — UGC 3995 
(Seyfert 2) and NGC 1167 (LINER) — do not exhibit 
detectable HI disturbances. In the optical, UGC 3995 
(Seyfert 2) forms an overlapping pair with UGC 3995A, 
and may therefore comprise an interacting system. 
NGC 1167 (LINER) is the only galaxy in our sample 
where we find no evidence whatsoever for disturbances 



in either the optical or HI. 

Markarian 993 

Markarian 993 (Mrk 993) is a Sa galaxy hosting a 
Seyfert 1.5 nucleus. No optical disturbances nor bars 
are visible in this galaxy even at the angular resolution 
of the HST (Malkan, Gorjian, & Tarn, 1998; Pogge & 
Martini, 2002). 

In our HI image (Fig. 24), we detected Mrk 993 
as well as a comparatively optically-dim galaxy lying 
— 8'.7 (—175 kpc) to the east, hereafter referred to as 
HI 0126+321. These galaxies have systemic velocities 
that differ by —450 km s" 1 , the largest of any confirmed 
or candidate interacting galaxy pairs in our sample. In 
the moment maps (Fig. 24), the HI disk of Mrk 993 ap- 
pears to be quite normal except for a weak protuberance 
on its northern-eastern side. This protuberance can also 
be seen in the channel maps (Fig. 25), and points east 
in the direction of HI 0126+321. Weak protuberences 
pointing west, in a direction opposite to HI 0126+321, 
also can be seen on the northern-eastern side of the HI 
disk. We tentatively associate these weak disturbances 
to tidal interactions with HI 0126+321. 

UGC 3157 

UGC 3157 is a barred Sbc galaxy hosting a Seyfert 2 
nucleus. It has a highly asymmetric appearence in the 
DSS2 image caused by a relatively dim and diffuse swirl 
on the eastern side of the galaxy. There are, however, no 
previously cataloged optical galaxies nor are there any 
obvious galaxies in the DSS2 image within several ar- 
cminutes of UGC 3157. 

Our HI image of UGC 3157 (Fig. 27) reveals that this 
galaxy has an asymmetric HI disk. The centroid in inte- 
grated HI intensity is shifted east of the geometrical and 
kinematic center of the HI disk, as well as the center of 
the optical galaxy. An examination of the channel maps 
(Fig. 28) reveal that the HI emission is stronger on the 
eastern side of the disk, tracing the diffuse optical swirl 
mentioned above. The disturbance seen in the optical 
is therefore also apparent in HI, albeit less clearly pre- 
sumably because of our lower angular resolution. The HI 
kinematics of the galaxy, however, appears normal. 

We also detected in HI a relatively dim object located 
-4f8 (-95 kpc) to the north-west of UGC 3157, differ- 
ing in systemic velocity by —150 km s _1 . In the DSS2 
image, this object appears to be diffuse, apart from a 
relatively bright star-like object at its north-western 
periphery, and hence is a genuine galaxy. It is henceforth 
referred to as HI 0446+185. The HI kinematic axis 
of this galaxy coincides with its optical major axis, 
providing additional confidence that we have correctly 
identified the optical counterpart of HI 0446+185. 
It is possible that the optical and HI disturbances 
seen in UGC 3157 is caused by tidal interactions with 
HI 0446+185. 

Markarian 1419 

Markarian 1419 (Mrk 1419) is a Sa galaxy hosting a 
LINER nucleus. It has an asymmetric appearence in 
the DSS2 images, caused by a relatively dim and diffuse 
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swirl on the north-eastern side of the galaxy. There are, 
however, no cataloged galaxies in the vicinity of Mrk 1419 
at a comparable rcdshift. 

Our HI image of Mrk 1419 (Fig. 29) reveals that this 
galaxy has an asymmetric HI disk. Analogous to the 
case of UGC 3157, in Mrk 1419 the centroid in integrated 
HI intensity is shifted north-east of the geometrical 
and kinematic center of the HI disk, as well as center 
of the optical galaxy. In the channel maps (Fig. 30), 
the HI emission is stronger and clearly disturbed on 
the north-eastern side of the disk, coincident with the 
diffuse optical swirl mentioned above. The disturbance 
seen in the optical is therefore also apparent in HI, 
albeit less clearly presumably because of our lower 
angular resolution. The HI kinematics of the galaxy, 
however, appears normal, just as in the case UGC 3157. 
Compared with UGC 3157, the only difference is that 
we did not detect in HI any neighboring galaxy (nor any 
other galaxies in the same field of view) that could have 
caused the disturbances seen in Mrk 1419. 

IRAS 14082+1347 

IRAS 14082+1347 has an uncertain Hubble type in 
Hyperleda. This galaxy hosts a LINER nucleus. There 
are no cataloged galaxies at a comparable redshift, nor 
are there any obvious galaxies in the DSS2 images, in the 
vicinity of IRAS 14082+1347. 

We detected IRAS 14082+1347 at a relatively low sig- 
nal to noise (Fig. 31). This galaxy has a HI gas mass of 
only ^3.5 x 10 7 M , which is at least an order of mag- 
nitude smaller compared with the other galaxies in our 
sample. The HI body of this galaxy is clearly displaced 
from the optical center of the galaxy. In addition, the HI 
major and apparently also kinematic axes of this galaxy 
is clearly different from its optical major axis. 

Hcislcr & Vader (1994) reported shells around the 
amorphous body of IRAS 14082+1347. Given its 
appearence in the DSS2 image, and the HI properties 
of this galaxy, it is possible that IRAS 14082+1347 is 
an early-type (elliptical) galaxy that has cannibalized a 
gas-rich dwarf galaxy. 

UGC 3995 

UGC 3995, a Sb galaxy, is one of an optically- 
overlapping pair of spiral galaxies separated by just ~0.'5 
(~11 kpc) lying at the same redshift. UGC 3995 is the 
brighter eastern galaxy, and hosts a Seyfert 2 nucleus. 
The western galaxy is UGC 3995A, a Sbc galaxy. (We 
take the names of the galaxies from the NED. In some lit- 
erature sources, UGC 3995 is referred to as UGC 3995A 
(as in Veron-Cetty & Veron, 1998), and UGC 3995A re- 
ferred to as UGC 3995B.) Despite their overlapping disks, 
Marziani et al. (1999) found no signs of interactions be- 
tween the two galaxies from their HST image. 

This pair of galaxies cannot be separated in our HI 
image (Fig. 33). An examination of the channel maps 
(Fig. 34) indicate that we have detected HI emission 
from both galaxies. In both cases, the HI traces nor- 
mal rotating disks with no obvious signs for disturbances. 

NGC 1167 



NGC 1167 is a SO galaxy that hosts a LINER nucleus. 
The galaxy does not appear to be disturbed in the opti- 
cal. 

Our HI image of NGC 1167 (Figs. 35 and 36) reveals 
an enormous gas disk around this galaxy, the largest seen 
in our sample. The HI kinematics are consistent with a 
normal rotating disk with no signs of disturbances. In 
addition to HI in emission, we also detected spatially- 
unresolved HI in absorption (spanning a velocity width 
of ^250 km s _1 ) towards the center of the galaxy. Apart 
from NGC 1167, we detected in HI a relatively dim 
galaxy henceforth referred to as HI 0302+352 lying ^6'9 
(147 kpc) to the north-cast of NGC 1167. This galaxy 
also is not obviously disturbed in the optical nor HI. 

4.2. Other Galaxies detected in Target Fields 

As described in § 4.1, we detected a number of galaxies 
in HI around our target objects, many of which we show 
to be definitely (in Group I), probably (in Group II), 
or possibly (in Group III) interacting with the target 
galaxies. In addition, we detected a number of other 
galaxies mostly at larger projected separations from the 
target galaxies than the (candidate) interacting neigh- 
boring galaxies. The optical and HI images (moment 
maps) of these other galaxies are shown in Figures 38- 
43. Many of these galaxies, like the (candidate) interact- 
ing neighboring galaxies, are relatively faint (and small) 
in the optical, and are not obvious candidates for neigh- 
boring galaxies. The properties of all the neighboring 
galaxies are listed in Tables 4 and 5. We separate, in 
our discussion, those neighboring galaxies lying within 
the ~80% response (radius of ~9f 1 or ^190 kpc) of the 
primary beam and hence observed at approximately the 
same sensitivity as the active galaxies, and those lying 
further out up to the ~50% (radius of ~15f8 or ^320 kpc) 
response of the primary beam and hence observed at a 
significantly lower sensitivity. 

Altogether, we detected nineteen surrounding galaxies 
within ~190 kpc of our ensemble sample. Ten of these are 
confirmed interacting neighboring galaxies in Group I, 
three are probable interacting neighboring galaxies in 
Group II, and two are possible interacting neighboring 
galaxy in Group III (i.e., totaling fifteen of the nineteen). 
Thus, the vast majority of galaxies detected in HI within 
~190 kpc of the target galaxies are confirmed or likely 
to be interacting with the target galaxies. An additional 
eight galaxies lie between ^190 kpc and ^320 kpc of the 
target galaxies, two of which (in the range 200-210 kpc) 
we show to be definitely (Group I) interacting with the 
target galaxies. 

In Figure 44 (upper panel), we show the number as well 
as cumulative fraction of Seyfert galaxies in our ensemble 
sample having (candidate) interacting neighboring galax- 
ies (i.e., in Groups I, II, or III), plotted as a function 
of their projected separations (for the few with multi- 
ple interacting neighboring galaxies, their nearest such 
neighbor). Within 210 kpc, there are only two neigh- 
boring galaxies not identified to be (possibly) interacting 
with the same Seyfert galaxy, comprising two of the three 
neighboring galaxies (the other being identified as the 
interacting neighbor) around ARK 539; including these 
galaxies do not therefore change the statistics. A dashed 
line is drawn at a projected separation of 190 kpc, within 
which our HI detection threshold is relatively uniform. 
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As can be seen, the fraction of Seyferts with neighboring 
galaxies increases steeply with projected separation up 
to ~100 kpc; by this separation, ~80% of the Seyferts 
already have neighboring galaxies. Within a projected 
separation of ~180 kpc, nearly all (~95%) have neigh- 
boring galaxies. 

In Figure 44 (lower panel), we plot the difference in 
radial velocity between the Seyfert and all their neigh- 
boring galaxies galaxies detected in HI within 210 kpc. 
We use different symbols for those in the different inter- 
acting classes (e.g., our defined Group I, II, and III), as 
well as those not apparently interacting with the Seyfert 
galaxy. The vast majority of the confirmed or probable 
interacting companion galaxies lie within ±100 km s _1 , 
and all but one within about ±200 km s _1 . 

4.3. Ensemble Statistics 

We now compile statistics relevant to this study from 
the results for the ensemble sample of twenty-three galax- 
ies that we mapped in HI. Two are Seyfert Is, three are 
Seyfert 1.5s, two are Seyfert 1.9s, nine are Seyfert 2s, 
two are unclassified Seyfert types (i.e., total of eighteen 
Seyferts in all), four are LINERs, and one is a H II galaxy. 

4.3.1. Comparison with optical starlight 

Based on DSS2 images, six of the twenty-three galaxies 
(~26%) — comprising five Seyferts and one LINER — 
are disturbed in the optical. Among these six, only three 
(NGC 7679, Mrk 1040, and Mrk 341, all Seyferts) ex- 
hibit optical tidal features that can be confidently (in the 
case of NGC 7679, with the aid of our HI image) traced 
to interactions with neighboring galaxies. Two of these 
galaxies (Mrk 1040 and Mrk 341) comprise the closest 
interacting pairs, and overlap in the optical with their in- 
teracting neighbour. Two (NGC 7679 and Mrk 1040) ex- 
hibit corresponding HI tidal features, and one (Mrk 341, 
itself not detected in HI) corresponding HI tidal features 
in its interacting companion. The remaining three galax- 
ies that appear to be disturbed in the optical (UGC 1395, 
UGC 3157, and Mrk 1419) all exhibit asymmetric disks 
in the form of an extended spiral arm. Two (UGC 1395 
and UGC 3157) are Seyferts, and one (Markarian 1419) 
is a LINER. All three also exhibit corresponding HI dis- 
turbances. 

The remaining seventeen of the twenty-three galaxies 
(~74%) are not visibly disturbed in DSS2 images. One 
(NGC 7469) lies in an interacting system where only its 
neighbour exhibits tidal features. Another (NGC 5548) 
cxibits an extended feature detectable only in a deep 
ground-based optical image that we confirm from our 
HI image to be a tidal feature. In dramatic contrast 
with their optical appcarences, fifteen of these seventeen 
galaxies are disturbed in HI. Ten (Mrk 1, NGC 7682, 
ARK 539, NGC 5548, Mrk 461, Mrk 1157, NGC 841, 
Mrk 573, NGC 7591, and NGC 7469), comprising nine 
Seyferfs and one LINER, exhibit HI tidal features that 
can be directly traced to interactions with neighbor- 
ing galaxies (i.e., those in Group I). The remaining five 
(MS 04595±0327, NGC 513, Mrk 1158, Mrk 993, and 
IRAS 14082±1347) all exhibit spatial and kinematic dis- 
turbances in HI, and which in four cases is likely pro- 
duced by interactions with neighboring galaxies detected 
also in HI. Three are Seyferts, one a Liner, and one a 
H II galaxy. Only two of the twenty-three galaxies in our 



ensemble sample do not exhibit any visible disturbaces 
in (optical and) HI. One is a Seyfert, and the other a 
LINER. 

In summary, only six of the twenty-three galaxies 
(^26%) in our ensemble sample are visibly disturbed in 
DSS2 images. In dramatic contrast, twenty-one of the 
same twenty-three galaxies (~91%) are disturbed in HI, 
thirteen (^62%) of which are interacting with neighbor- 
ing galaxies (Group I), another four (~19%) probably in- 
teracting with neighboring galaxies (Group II), and the 
remaining four (^19%) weakly disturbed (Group III). 
In nearly all these cases, the HI is disturbed not just 
spatially but also kincmatically on global scales. Of 
the eighteen galaxies in our ensemble sample classified 
as Seyferts, only five (~28%) are visibly disturbed in 
DSS2 images. By contrast, seventeen (~94%) exhibit 
HI disturbances, twelve (^67%) of which are interact- 
ing with neighboring galaxies (Group I), another three 
(~17%) probably interacting with neighboring galaxies 
(Group II), and the remaining two (^12%) weakly dis- 
turbed (Group III). 

Our parent sample contains all the Seyfert galaxies in 
the CfA Seyfert and Markarian catalogs within our se- 
lected redshift range (see §2). The trends seen in our 
ensemble sample are also seen in the complete subsam- 
ples from these catalogs. Of the seven Seyfert galaxies 
in our sample present also in the CfA Seyfert catalog 
(Huchra ct al., 1983), only one (~14%) of these galaxies 
is visibly disturbed in DSS2 images. By contrast, in HI, 
five (~71%) are found to be interacting with neighbor- 
ing galaxies (Group I), and the other two (^29%) dis- 
turbed (Group III) . Fifteen of the galaxies in our sample 
are present in the Markarian catalog (Markarian et al., 
1989), twelve of which (including two classified here as 
LINER and H II respectively) are in our ensemble sam- 
ple. Four of these twelve galaxies (~33%) are visibly 
disturbed in DSS2 images, whereas nine (~75%) are in- 
teracting (Group 1) and three (~25%) probably or weak 
disturbed (Group II or III) in HI. Among just the ten 
galaxies classified here as Seyferts, three (~30%) are vis- 
ibly disturbed in DSS2 images, whereas nine (~90%) are 
interacting with neighboring galaxies (Group 1) and one 
(-10%) weakly disturbed (Group III) in HI. 

4.3.2. Comparison with other galaxies in target fields 

The other galaxies detected in our target fields, lying 
within a projected separation of ~320 kpc from the tar- 
get objects, may provide a natural sample that is best 
matched in mean environment to our ensemble sam- 
ple. These galaxies, however, are clearly not matched 
in optical luminosities. Where cataloged and have ap- 
parent _B-band magnitudes listed in Hyperleda, many 
have corresponding absolute _B-band magnitudes close 
to or spanning the range between the Small and Large 
Magellanic Clouds. (If placed at the same distance, the 
SMC would have a corrected apparent S-band magni- 
tude ttib ~ 18.5, and the LMC tub ~ 16.5.) Those not 
cataloged or do not have apparent -B-band magnitudes 
listed in Hyperleda likely have even smaller optical lu- 
minosities. These galaxies also are not matched in other 
optical properties such as size or Hubble type; where cat- 
aloged in Hyperleda, many are later-type Sb-c galaxies. 

Eleven of the nineteen neighboring galaxies (^58%) 
detected at approximately the same sensitivity in HI as 
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the ensemble sample (i.e., those lying within ~190 kpc 
of the target galaxies) are disturbed in HI. All eleven 
of these disturbed neighboring galaxies are among the 
fifteen we identify to be definitely (Group I), probably 
(Group II), or possibly (Group III) interacting with the 
target galaxies. The trend therefore is towards a higher 
incidence of HI disturbances among the target compared 
with their neighboring galaxies. We caution, however, 
that this may be an observational selection effect, as the 
HI gas masses and extents of the neighboring galaxies are 
usually much smaller than that of the target galaxies, and 
hence any comparable level of disturbances more difficult 
to detect. 

4.3.3. Tidal interactions versus mergers 

Our HI images allow us to directly determine the cause 
of the gaseous disturbances in many of the target galax- 
ies. Thirteen of the galaxies in Group I (twelve Seyferts 
and one LINER) out of the ensemble sample of twenty- 
three galaxies (~57%) have recently experienced tidal 
interactions with neighboring galaxies. Adding to this 
the three galaxies in Group II and the two galaxies in 
Group III (all Seyferts) that have likely been disturbed 
by neighbouring galaxies also detected in HI, eighteen of 
the twenty-three galaxies in our ensemble sample (~78%) 
have likely experienced recent tidal interactions with 
neighboring galaxies. Among those classifed as Seyferts, 
twelve of eighteen (~67%) have experienced recent tidal 
interactions with neighboring galaxies (Group I), and an- 
other five (~28%) likely or probably have experienced 
such interactions (Group II or III); i.e., in total, seven- 
teen of the eighteen (~94%) Seyfert galaxies. 

Only three of the twenty-one galaxies that are dis- 
turbed in HI do not have neighboring galaxies cataloged 
in the optical or detected in HI. One is a H II galaxy, 
and the remaining two LINERS. Thus, only three of 
the twenty-one galaxies in our ensemble sample (~14%) 
are good candidates for recent (minor) mergers, none of 
which are Seyferts. Note that another two galaxies (one 
Seyfert and one LINER) are not disturbed in either the 
optical or HI (and for which we cannot immediately rule 
out the possibility of ancient minor mergers). Of course, 
in all these cases, the disturbances or interacting neigh- 
boring galaxies may simply not have been detectable in 
our HI observations. 

5. INTERPRETATION 

5.1. Bias towards Interacting Systems? 

Given the prevalence of tidal interactions or distur- 
bances seen in our HI images, is it possible that our sam- 
ple is somehow biased towards interacting systems? This 
could possibly happen if we had imposed in our selection 
criteria a threshold in HI gas mass based on single-dish 
observations, thus ensuring that all the galaxies are de- 
tected. The more gas-rich systems could preferentially 
comprise interacting or merging systems; i.e., contain two 
or more galaxies in the single-dish beam, or constitute 
merger remnants. We remind the reader, however, that 
our selection is entirely blind to the HI content of the 
galaxies; indeed, the HI gas mass of the detected Seyfert 
galaxies span the broad range ~10 8 -10 10 M©, with a me- 
dian of ^6 x 10 9 M©. For the same reason, if only a 
relatively small number of the more gas-rich galaxies are 



detected, this could skew the results towards interacting 
systems. As mentioned at the start of §4, only two of 
our sample galaxies were not detected (with two others 
omitted because of RFI) ; in one case (Mrk 359) the pres- 
ence of a strong continuum source resulted in a higher 
noise detection threshold than targeted, and in the other 
case (2E 0507+1626) the redshift is not accurately known 
(see Greene, Lim, & Ho, 2004). Because we detected the 
vast majority of the galaxies in our sample, our results 
arc unlikely to be biased towards the more gas-rich and 
therefore perhaps preferentially interacting systems. 

Is there any evidence from their optical morphology 
that we had somehow, perhaps just because of small 
number statistics, preferentially selected interacting sys- 
tems? Only five of the eighteen galaxies (~28%) in our 
ensemble sample classified as Seyferts are visibly dis- 
turbed in DSS2 images. Including the remaining four 
galaxies classified as LINERS and one galaxy classified as 
H II, six of the twenty-three galaxies (~26%) are visibly 
disturbed in DSS2 images. The relatively small fraction 
of optically disturbed — let alone (apparently) optically 
interacting — galaxies in our sample is typical of that 
found in other optical studies utilizing different samples 
and images. 

For instance, de Robertis, Hayhoe & Yee (1998) have 
studied thirty-four Seyfert galaxies selected from the CfA 
Seyfert catalog (nearly half of which can also be found 
in the Markarian catalog) using much deeper optical im- 
ages than DSS2 images. Only about 26% (9/34) of the 
Seyferts in their sample were classified as disturbed (a 
fraction comparable with that found in our sample) . The 
majority (7/9) of these disturbed Seyferts and their in- 
teracting neighbors form overlapping pairs, with some 
perhaps comprising late mergers. Most (28/34) of the 
Seyferts in that sample have absolute B-band magni- 
tudes (as listed in Veron-Cetty & Veron, 1998) spanning 
the same range as our sample. Considering only this sub- 
sample, about 29% (8/28) of the Seyferts are disturbed. 
These galaxies, however, occupy a much broader redshift 
range than our sample, spanning 0.0033 < z < 0.0422 
with an average redshift of z = 0.0189; i.e., reaching 
both lower and higher rcdshifts than our sample. Dis- 
turbances are of course more difficult to detect in more 
distant galaxies, although de Robertis, Hayhoe & Yee 
(1998) found no apparent distance bias in their ability to 
detect disturbances with the disturbed galaxies spread 
quite evenly over the entire redshift interval of their sam- 
ple. 

Although a direct comparison between the incidence of 
optical disturbances we and de Robertis, Hayhoe & Yee 
(1998) find in our respective samples cannot be made, 
we note that the disturbances and/or tidal features de 
Robertis, Hayhoe & Yee (1998) sec in all their disturbed 
galaxies can also be seen in DSS2 images. Thus, a study 
using just DSS2 images would reach the same conclusion 
as de Robertis, Hayhoe & Yee (1998) on the incidence 
of optical disturbances in their sample. There are only 
two galaxies in our sample in common with that of de 
Robertis, Hayhoe & Yee (1998), Mrk 461 and NGC 5548. 
de Robertis, Hayhoe & Yee (1998) find neither galaxies 
to be optically disturbed, whereas we find both to be 
interacting with neighboring galaxes detected also in HI 
(i.e., Group I). 
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5.2. Bias towards Higher Density Environments? 

We also checked to see whether our sample is, perhaps 
just by chance, somehow biased towards higher density 
environments compared with the general Seyfert popu- 
lation. Eleven of the twenty-three Seyfert galaxies in 
our parent sample were studied by Dahari (1984), about 
half of whose sample comprise Markarian galaxies. This 
study remains to date that which has found the largest 
disparity between the fraction of Seyfert galaxies likely 
to have physical companions (~15%) compared with a 
control sample of inactive galaxies (~3%). Only two of 
the eleven (~18%) Seyfert galaxies in our parent sample 
were listed by Dahari (1984) to have projected neighbor- 
ing galaxies, similar to that found in their overall sample 
of Seyfert galaxies. Nine (including the two with pro- 
jected neighboring galaxies) of these eleven galaxies are 
in our ensemble sample, all of which we find in our HI im- 
ages to be definitely or probably interacting with neigh- 
boring galaxies; i.e., seven (<~78%) are in Group I, and 
two (~22%) in Group II. 

Why did Dahari (1984) not find many of the interact- 
ing neighboring galaxies that we detected in HI? The lim- 
iting magnitude used in his work was 18.3, comparable 
with that of the SMC at the distance of our sample and 
hence sufficient to detect most of the neighboring galax- 
ies that we detected. Dahari (1984), however, searched 
for projected neighboring galaxies only at angular sep- 
arations within three times the optical diameters of the 
Seyfert (or control inactive) galaxies (the search radius 
therefore differs for Seyfert/inactive galaxies with dif- 
ferent optical diameters). Projected neighboring galax- 
ies (genuine neighboring galaxies, as it turns out) were 
found for NGC 513 and NGC 7469, but not for Mrk 1, 
352, 359, 573, 1157, NGC 5548, 7682, ARK 539, and 
UGC 1395. In all the latter cases (excluding Mrk 352 and 
Mrk 359, which are not in our ensemble sample) apart 
from UGC 1395, the actual interacting companion galaxy 
is located beyond the search radius used; the possible in- 
teracting companion galaxy of UGC 1395 (Group II) is 
unrecognizable even in modern (DSS2) optical images as 
a galaxy. 

Rafanelli et al. (1995) have searched for an ex- 
cess of projected neighboring galaxies around all (261) 
spectroscopically-classified Seyfert galaxies listed in the 
Veron-Cetty & Veron (1991) catalog at z < 0.05 (mostly 
at significantly lower redshifts than the upper limit) and 
6 > —23°. No significant differences were found between 
the frequency of Seyfert 1 and 2 galaxies with projected 
neighboring galaxies, but this frequency was found to be 
much higher than in matched samples of inactive field 
spiral galaxies. Twelve of the Seyfert galaxies in our par- 
ent sample are included in that study. Only four of these 
twelve (33%) Seyfert galaxies are listed in Rafanelli et al. 
(1995) as having projected neighboring galaxies, a frac- 
tion comparable with that found for their entire sample 
(~26%). Ten of these twelve galaxies arc in our ensemble 
sample, seven of which (70%) are interacting with neigh- 
boring galaxies detected also in HI (Group I) , two others 
(20%) probably interacting with neighboring galaxies de- 
tected also in HI (Group II), and the remaining galaxy 
(10%) probably disturbed in HI (Group III). 

Rafanelli et al. (1995) used the same search criteria 
as Dahari (1984), but required the Seyfert and pro- 



jected neighboring galaxy to have a difference in ap- 
parent magnitude not exceeding 3. They identified pro- 
jected neighboring galaxies for Mrk 1, 1040, NGC 7469, 
IRAS 1408+1347 (almost certainly the genuine neighbor- 
ing galaxies for the first three galaxies; we do not detect 
an interacting neighboring galaxy for IRAS 1408+1347, 
classified here as a LINER), but not for Mrk 352, 359, 
573, 993, 1157, NGC 513, 5548, and 7682. Once again, 
in the latter case (apart from NGC 513) the interact- 
ing neighboring galaxies detected in HI lie beyond the 
search radius used. Rafanelli et al. (1995), unlike Dahari 
(1984), failed to identify a projected neighboring galaxy 
for NGC 513; its interacting neighboring galaxy falls out- 
side the magnitude difference criteria used by Rafanelli 
et al. (1995). 

There is therefore no independent evidence from op- 
tical studies that our Seyfert sample is somehow biased 
towards higher density environments than the general 
Seyfert population. On the contrary, when we can find 
a large number of the Seyfert galaxies in our sample 
included in optical studies, the fraction with projected 
neighboring galaxies as identified in those studies is iden- 
tical with the fraction found in their overall sample. 

5.3. Prevalence of Tidal Interactions in HI Gas 

The results from our HI study, by contrast with those 
from optical studies, dramatically alter our perception of 
the frequency (and strength) of large-scale disturbances 
in local Seyfert galaxies. Of the eighteen galaxies in 
our ensemble sample of twenty-three galaxies classified 
as Seyferts, seventeen (~94%) exhibit HI disturbances, 
compared with only five (~28%) exhibiting optical dis- 
turbances in DSS2 images. In twelve (including Mrk 341) 
of these eighteen cases (67%), the HI disturbances can 
be directly traced to tidal interations with neighboring 
galaxies detected also in HI. In the other five of six cases, 
the observed HI disturbances are likely produced by tidal 
interactions with neighboring galaxies detected also in 
HI. Thus, nearly all (~94%) of the Seyfert galaxies in 
our ensemble sample can be directly or indirectly shown 
to be involved in recent or ongoing tidal interactions with 
neighboring galaxies. 

Our parent sample comprises all known (circa 1998) 
Seyfert (disk) galaxies within a given redshift interval in 
the northern hemisphere. These galaxies were compiled 
from the same surveys that found a large fraction of the 
known population, and are therefore likely representa- 
tive of this heterogeneously selected general population. 
The only intended selection bias is towards Seyfert nu- 
clei with relatively high optical luminosities; specifically, 
within about two orders of magnitude of the lower limit 
for QSOs. Our results therefore suggest that a large frac- 
tion (if not the vast majority) of local Seyfert galaxies 
with high nuclear luminosities are involved in or have 
recently experienced tidal interactions with neighboring 
galaxies. Much more often than not, these interactions 
are visible only in HI gas. The prevalence of such tidal 
interactions constitutes a shift in our understanding of 
local Seyfert galaxies, and is the first of the two impor- 
tant results of our study. 

We note that there are only two Seyfert 1 and three 
Seyfert 1.5, compared with eleven Seyfert 2, in our 
ensemble sample (the remaining two being unclassified 
Seyfert types). Given that many optical studies find an 



1G 



Kuo et al. 



excess of projected neighboring galaxies around Seyfert 2 
but not Seyfert 1 galaxies, it may be that these two 
classes of galaxies have genuinely different environments 
and that our results can be confidently extrapolated only 
to Seyfert 2 galaxies. On the other hand, all the Seyfert 1 
and Seyfert 1.5 galaxies in our sample also are interacting 
or disturbed in HI (three classified in our Group I, one in 
Group II, and one in Group III), and so (with the caveat 
of small number statistics) we do not see a dichotomy in 
the frequency of tidal interactions between Seyfert 1 and 
Seyfert 2 galaxies in HI. 

5.4. Implicating Tidal Interactions in Triggering AGNs 

The observed prevalence of tidal interactions in Seyfert 
galaxies does not, by itself, necessarily imply that such 
interactions are responsible for initating events that lead 
to their AGN activity. It may be that the HI disks of 
many galaxies — both active and inactive — are mod- 
erately to strongly disturbed at their outskirts by tidal 
interactions with neighboring galaxies. To test whether 
this is the case, we have conducted a HI imaging survey of 
a matched control sample of inactive galaxies at an iden- 
tical spatial resolution and detection threshold to those 
attained here. The results of that study are reported in 
a companion paper by Tang et al. (2007). Below, we 
provide just a synopsis of that work. 

Tang et al. (2007) selected a sample of twenty-seven 
inactive galaxies from the CfA Redshift Survey (Huchra 
et al., 1983, 1995) that were closely matched on a one-to- 
one basis in Hubble type, to the degree possible galaxy 
optical luminosity, and in range of size and inclination 
to our parent active galaxy sample. Other than that, 
the inactive galaxies were essentially selected at random 
from the field, and in this way we hoped representative 
of the environment of inactive field galaxies. To reduce 
the required observing time, Tang et al. (2007) selected 
inactive galaxies at redshifts of 0.0075 < z < 0.0085, 
half that of our active galaxy sample. All the galaxies 
in the control sample were imaged in HI at the same 
spatial resolution (~20 kpc) and detection threshold as 
our active galaxy sample. 

Twenty-one of the twenty-seven inactive galaxies de- 
tected in HI comprise the ensemble control sample. 
These galaxies have approximately the same distribution 
in HI gas masses as our ensemble sample of active galax- 
ies. Tang et al. (2007) found that only four of these 
twenty-one galaxies (^19%) exhibit spatial and usually 
also kinematic disturbances on galactic (> 20 kpc) scales. 
One of these disturbed galaxies has since been found to 
be a Seyfert 2 galaxy, and so only three of the twenty 
(15%) inactive galaxies in our ensemble control sample 
are actually disturbed in HI. By contrast, seventeen of 
the twenty-one (~81%) galaxies show no HI disturbances 
whatsover on galactic (> 20 kpc) scales. Excluding again 
the one disturbed galaxy since found to be a Seyfert 2 
galaxy, seventeen of the twenty (85%) inactive galaxies 
in our ensemble control sample are not disturbed in HI. 

The contrast between the incidence of HI disturbances 
in the active (or just Seyfert) and inactive galaxy sam- 
ples is dramatic, and difficult to explain away as selection 
or other spurious effects. In Tang et al. (2007), we com- 
pare in greater detail the contrasting results found for 
the Seyfert and control samples in our HI study with 
the results found for such samples in optical studies, 



and show that all these results need not be incompat- 
ible given the different criteria used to find candidate 
neighboring galaxies in optical studies. The prevalence 
of HI disturbances in our sample of active galaxies, but 
the low incidence of disturbances at a comparable level 
in our control sample of inactive galaxies, directly impli- 
cates tidal interactions in initiating events that lead to 
optically-luminous Seyfert activity in a large fraction of 
local disk galaxies. This constitutes another shift in our 
understanding of local Seyfert galaxies, and is the second 
of the two important results of our study. As before, our 
results can be more confidently extrapolated to Seyfert 2 
galaxies, although we find no difference in our results for 
Seyfert 1 galaxies. 

6. DISCUSSION 

If our sample is indeed representative of the general 
population, our results have interesting and important 
implications for the nature of gravitational interactions 
involving Seyfert galaxies, the timescale for gas to be 
transported inwards to fuel their central SMBHs, and 
the lifetime of optically-luminous Seyfert activity. 

6.1. Bound Interactions or Flybys? 

Do the observed tidal interactions constitute gravita- 
tionally bound interactions or flybys? Seyfert galaxies 
preferentially comprise luminous early-type disk galax- 
ies, as is the case also in our sample as can be seen in 
Table 1. To esimate whether the Seyfert galaxies in our 
ensemble sample arc gravitationally bound to their con- 
firmed or likely interacting neighboring galaxies, let us 
assume that these Seyfert galaxies have a typical mass 
(including dark matter) of ~10 12 M , comparable with 
that of our own Milky Way galaxy. The projected spa- 
tial separation between the Seyferts and their (candi- 
date) interacting neighbors ranges from a minimum of 
~6 kpc to a maximum of ~210 kpc, with the majority 
(~85%) lying within ^100 kpc. The projected escape ve- 
locity (i.e., velocity along the plane of the sky) therefore 
ranges from ^1200 km s -1 at the closest separation to 
^200 km s _1 at the largest separation, and corresponds 
to ^800 km s _1 at ^100 kpc. For comparison, most of 
the interacting neighboring galaxies differ in radial veloc- 
ities by < 100 km s _1 , and nearly all by < 200 km s _1 , 
from their respective Seyfert galaxies. Scaling as the 
square root of the Seyfert galaxy mass, the escape ve- 
locity remains in most cases larger than the observed 
difference in radial velocities even if the Seyfert galaxy 
mass is an order of magnitude smaller than has been as- 
sumed. Because the overall difference in radial velocities 
should be statistically comparable with the overall differ- 
ence in velocities along the plane of the sky, most of the 
interacting neighboring galaxies arc likely to be gravita- 
tionally bound to their respective Seyfert galaxies. 

6.2. Gas Infall Timescale 

The relative simplicity of many of the tidal features 
observed — clearly identifiable tidal bridges and tails — 
suggests that most of the systems studied are caught at 
a relatively early phase of an encounter; i.e., seen just 
after their first close approach. If this is indeed the 
case, then our results place a stringent upper limit on the 
timescale for gas to be transported inwards to fuel the 
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central SMBH. For our ensemble sample of Seyfert galax- 
ies at a mean redshift of z — 0.016 (distance —71 Mpc), 
the time, t, since closest approach (when tidal forces are 
strongest, and during which theoretical simulations sug- 
gest tidal distortions first become appreciable) in an in- 
teracting pair is: 

where ^pc is the spatial separation (in kpc) between the 
two galaxies, and v is their relative velocity in the plane 
of the sky. As mentioned above, the projected spatial 
separation between the (candidate) interacting neighbor- 
ing galaxies and the Seyfert galaxies in our ensemble sam- 
ple ranges from a minimum of —6 kpc to a maximum of 
-210 kpc, with most (-85%) within -100 kpc. If we 
assume that their relative velocities along the plane of 
the sky is comparable with their relative radial veloci- 
ties of usually < 100 km s _1 , then the time since closest 
approach is — 10 8 -10 9 yrs. This is comparable with the 
orbital period in the outer regions of normal disk galax- 
ies, allowing perhaps sufficient time for inf ailing gas even 
at quite large radial distances from center to reach the 
nucleus. 

Numerical simulations show that mergers between two 
gas-rich disk galaxies with comparably large masses (i.e., 
major mergers) can be highly effective at transporting 
large amounts of gas into the centers of the merger rem- 
nants (e.g., Barnes & Hernquist, 1992, 1996; Mihos & 
Hcrnquist, 1996b). The same also can be true for minor 
mergers, where a relatively massive galaxy cannibalizes a 
lower mass but gas-rich galaxy (e.g., Mihos & Hcrnquist, 
1996a). On the other hand, while tidal interactions can 
strongly perturb the outskirts of galaxies, it is not imme- 
diately obvious how such interactions promote gas infall 
in the inner regions of the galaxy. More theoretical sim- 
ulations of such tidal interactions are needed to address 
the possible gas infall mechanisms involved. 

6.3. Longevity of Seyfert Activity 

If indeed most systems have just experienced their first 
close encounter, then the inferred timcscale since clos- 
est approach of 10 8 -10 9 yrs also places a stringent up- 
per limit on the longevity of Seyfert activity produced 
during a given interaction episode. If many Seyfert nu- 
clei shine continuously or sporadically for much longer 
timcscalcs, then we should sec a high fraction of sys- 
tems with relatively complex tidal features or in (near) 
mergers. Instead, there are very few good candidates 
for mergers, with most of the interacting systems hav- 
ing well separated neighboring galaxies and quite simple 
tidal features. 

7. SUMMARY AND CONCLUSIONS 

We reported in this paper the first systematic imag- 
ing program in atomic hydrogen (HI) gas of a represen- 
tative sample of Seyfert galaxies. Our parent sample 
comprised all twenty-seven disk galaxies in the Veron- 
Cetty & Veron (1998) catalog, plus another in the Vcron- 
Cetty & Veron (2000) catalog, classified as either Seyfert, 
LINER, or H II galaxies (the LINER and H II some- 
times classified as Seyferts in other catalogs or studies) 
at redshifts 0.015 < z < 0.017 (i.e., essentially uni- 
form distances) and with absolute -B-band magnitudes 



—19 > Mb > —23 (i.e., within about two orders of mag- 
nitude of the minimum optical luminosities for QSOs, 
as confirmed for over half with measured [OIII] lumi- 
nosities). We successfully imaged twenty-three of these 
galaxies in HI at essentially uniform sensitivities and spa- 
tial resolutions (—20 kpc), comprising eighteen Seyferts, 
four LINERs, and one H II galaxy. These twenty-three 
galaxies comprised our ensemble sample from which we 
drew the folio wing statistical results. 

1. In DSS2 images, only —26% (6/23) of the galaxies 
are visibly disturbed in the optical. Among those 
classified as Seyferts, only —28% (5/18) are vis- 
ibly disturbed in the optical. These results are 
comparable with other studies involving different 
Seyfert samples and using different optical images, 
all of which show that only a small minority of such 
galaxies are visibly disturbed in optical starlight. 

2. In dramatic contrast with their appearences in op- 
tical starlight, —91% (21/23) of the galaxies are dis- 
turbed in HI gas. The HI gas is disturbed not just 
spatially, but usually also kinematically, on global 
(i.e., galactic) scales. Among those classified as 
Seyferts, -94% (17/18) exhibit HI disturbances. 

3. The disturbances in —57% (13/23) of the galax- 
ies can be directly traced to tidal interactions with 
neighboring galaxies detected also in HI. Among 
those classified as Seyferts, the fraction is —67% 
(12/18). In another five galaxies (all Seyferts), the 
observed HI disturbances are most likely caused 
by tidal interactions with neighboring galaxies de- 
tected also in HI. Including the latter, —78% 
(18/23) of the galaxies in our ensemble sample, and 
—94% (17/18) of those classified as Seyferts, can be 
directly or indirectly shown to have recently experi- 
enced tidal interactions with neighboring galaxies. 

4. The confirmed or probable interacting neighbor- 
ing galaxies typically have projected separations 
of < 100 kpc and radial velocities differing by 
< 100 km s _1 from their respective Seyfert galax- 
ies. Many have optical B-band luminosities rang- 
ing from that of the Small to Large Magellanic 
Clouds. 

5. Three galaxies that are disturbed in HI (two LIN- 
ERs and one H II galaxy) do not have neighbor- 
ing galaxies apparent in the optical or detected in 
HI. These galaxies are the only good candidates for 
(minor) mergers. 

6. Only —9% (2/23) of the galaxies exhibit no de- 
tectable disturbances in HI, nor in the optical. One 
is a Seyfert, and the other a LINER. The Seyfert 
galaxy with no HI disturbances forms an optically- 
overlapping pair with a disk galaxy detected also 
in HI. 

Our results, if found to prevail in studies incorporating 
larger samples, have three important implications. 

1. The majority of local Seyfert galaxies with rela- 
tively high nuclear luminosities are involved in tidal 
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interactions with gravitationally bound neighbor- 
ing galaxies. Much more often than not, the re- 
sulting tidal features are visible in HI gas but not 
in optical starlight (in DSS2 images). 

2. The incidence of HI disturbances in a matched con- 
trol sample of inactive galaxies is only <~15% (Tang 
et al. , 2007) . The prevalence of HI disturbances in 
our Seyfert sample but low incidence of comparable 
HI disturbances in the control sample implicates 
tidal interactions in initiating events that lead to 
optically-luminous Seyfert activity in a large frac- 
tion of local field disk galaxies. 

3. The relative simplicity of the tidal features ob- 
served suggest that many Seyfert galaxies are in an 
early phase of an encounter with neighboring galax- 
ies; i.e., seen just after the first close approach. 
If this is indeed the case, then both the timescale 
for gas infall to fuel the central supermassive black 
hole, and the lifetime of Seyfert nuclei in their op- 
tically luminous phase, is < 10 8 10 9 yrs. 

Our work demonstrates in a systematic way the power 
of HI imaging observations to reveal tidal interactions 
where none are otherwise visible, and in most cases to 
positively identify and in many of the remaining cases 
pick out the most likely interacting neighboring galaxies. 
The application of this technique to Seyfert galaxies has 



provided the first strong link between tidal interactions 
and luminous Seyfert activity. 
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Note. — Units ol right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutcs, 
and arcscconds. Spectral type S denotes unclassified Seyfert galaxy. All data apart from Hubble types and inclinations 
are from the Veron-Cetty & Veron (1998) catalog, except for UGC 3157 from the Veron-Cetty & Veron (2000) catalog. 
UGC 3995 is referred to in Veron-Cetty & Veron (1998) as UGC 3995A. Hubble types and inclinations are from 
the Hyperleda database, where E? and S? denote unclassified elliptical and spiral galaxies respectively. The [OIII] 
luminosities are computed from measurements by Dahari & de Robertis (1988) and/or Whittle (1992) assuming H = 
50 km s~ x Mpc -1 and qo = 0. In the last column, reference 1 refers to Dahari & dc Robertis (1988) and refrence 2 to 
Whittle (1992). 
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TABLE 2 
Complete Subsamples 



Name 
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Note. — Those galaxies in 
our sample included in the CfA 
Seyfert catalog (Huchra ct al., 
1983) and/or Markarian cata- 
log (Markarian et al., 1989), la- 
belled in column 2 as CfA and 
Mrk respectively. 
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Observing Parameters 







Phase Center 




Flux 


Secondary 


Bandpass 


Observing 


Name 




R.A. 




Decl. 


calibrator 


calibrator 


calibrator 


date 


K4cirkaria.il 34 1 


00 


37 


16 


.6 


24 


02 


16 


01 37+331 


0029+349 


0029+349 


2000/7/28- 8/20 


Markarian 352 


00 


59 


53 


.3 


31 


49 


37 


0137+331 


0119+321 


0119+321 


2004/7/2 


Markarian 1 


01 


16 


04 


.7 


33 


09 


53 


0137+331 


0137+331 


0137+331 


2000/8/21 


NGC 513 


01 


24 


26 


.8 


33 


47 


58 


0137+331 


0137+331 


0137+331 


2001/10/30 


Markarian 993 


01 


25 


34 


5 


32 


04 


24 


0137+331 


0137+331 


0137+331 


2000/8/20 


Markarian 359 


01 


27 


19 


.5 


19 


09 


35 


0137+331 


0129+236 


0129+236 


2004/7/2 


M^arkarian 1157 


01 


33 


20 


.6 


35 


41 


42 


0137+331 


0137+331 


0137+331 


2001/10/30 


Markarian 1158 


01 


34 


59 


5 


34 


58 


46 


0137+331 


0137+331 


0137+331 


2001/11/2 


Markarian 573 


01 


43 


57. 


,8 


02 


20 


60 


0137+331 


0149+059 


01 37+331 


2004/7/2 


UGC1395 


01 


55 


24 


.5 


06 


43 


56 


0137+331 


0141+138 


0141+138 


2000/8/20 


NGC 841 


02 


11 


17 


.3 


37 


29 


49 


01 37+331 


0336+323 


01 37+331 


2000/8/20 


Markarian 1040 


02 


28 


35 





31 


18 


60 


0137+331 


0137+331 


0137+331 


2000/7/28; 8/20 


NGC 1 1 67 


03 


01 


42. 


.4 


35 


12 


21 


01 37+331 


0336+323 


0336+323 


2000/8/20 


IJGC3157 


04 


46 


48. 


9 


18 


24 


53 


01 37+331 


0521+166 


0521+166 


2000/8/21 


MS 04595+0327 


05 


02 


09 





03 


31 


51 


0137+331 


0503+020 


0503+020 


2001/10/30 


2E 0507+1626 


05 


10 


27 


.9 


16 


28 


22 


0137+331 


0521+166 


0521+166 


2001/11/2 


UGC3995 


07 


44 


28. 


.4 


29 


17 


06 


0137+331 


0741+312 


0137+331 


2000/8/20; 8/21 


Markarian 1419 


09 


40 


36 


A 


03 


34 


37 


1331+305 


1008+075 


1331+305 


2003/2/24 


Markarian 461 


13 


47 


17 


.7 


34 


08 


55 


1331+305 


1331+305 


1331+305 


2003/2/24 


IRAS 14082+1347 


14 


10 


41 


.4 


13 


33 


28 


1331+305 


1347+122 


1331+305 


2003/2/21 


Markarian 1510 


14 


17 


58. 


,6 


26 


24 


47 


1331+305 


1331+305 


1331+305 


2004/7/8 


NGC 5548 


14 


17 


59. 


5 


25 


08 


12 


1331+305 


1331+305 


1331+305 


2003/2/19 


ARK 539 


18 


28 


52 


.0 


50 


28 


48 


0137+331 


1829+487 


1829+487 


2000/8/20; 8/21 


NGC 7469 


23 


03 


27 


.7 


08 


54 


14 


0137+331 


2312+093 


2312+093 


2001/10/30 


NGC 7591 


23 


17 


58. 


6 


06 


33 


60 


0137+331 


2255+132 


2255+132 


2000/7/28; 8/21 


NGC 7679 


23 


28 


43 


.9 


03 


27 


02 


0137+331 


2255+132 


2255+132 


2000/7/28; 8/21 


NGC 7682 


23 


28 


43 


.9 


03 


27 


02 


0137+331 


2255+132 


2255+132 


2000/7/28; 8/21 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are 
degrees, arcminutes, and arcseconds. The secondary calibrators, located close to the individual target 
sources, are used for amplitude and phase calibration. 
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TABLE 4 

HI Gas Properties of Active Galaxies in Ensemble Sample 
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III 
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817 ±27 


4817 
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II 
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UGC 3995 


6170+367 
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6.67+0.39 


no 


III 


no 




Markarian 1419 


1681+271 
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1.81 ±0.27 


weak 


III 


yes 


Mrk 


Markarian 461 


1878+134 


4850 


2.31+0.16 


yes 
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III 
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NGC 5548 
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I 
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ARK 539 
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5123 
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NGC 7469 


8700+38 


4923 


10±1.2 
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NGC 7591 


17240+403 


4919 


23.98±0.55 


yes 
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6366+276 


5129 


8.85±0.34 


yes 
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yes 




NGC 7682 


7230+328 


5066 


10.05±0.46 


yes 


I 
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Cfa 



Note. — Assignment of groups based on observed spatial-kinematic distribution of HI gas (see text): Group I 
are those that show HI tidal features tracing interactions with neighboring galaxies (which also are detected in 
HI); Group II are those that show both spatial and kinematic disturbances in HI, but for which the observed 
HI disturbances cannot be directly (but can often be indirectly) traced to tidal interactions with neighboring 
galaxies; and Group III are those that show either weak or no HI disturbances. Optical disturbances assessed 
based on DSS2 images, and for some confirmed from HI maps. In the last column, those galaxies remarked with 
"Mrk" or "Cfa" comprise a complete subsamplc of Seyfert galaxies in either Markarian (Markarian et al., 1989) 
or CFA catalogue (Huchra et al., 1983) within out selected redshift interval. Note that Mrk 352, Mrk 359, and 
Mrk 1510 are not included in this table because these data were discarded due to cither nondetection or RFI. 
The exclusion of these three galaxies makes the Markarian subsamplc less "complete". 
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TABLE 5 

Other. Detected Galaxies in Active Galaxy Field 
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Sc 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, 
arcminutcs, and arcseconds. Coordinates are from the Hyperleda database except for UGC 3995A (from 
NED) and those not cataloged in the optical (from our HI maps). Optical rcdshifts are from NED where 
available. Apparent corrected B-band magnitudes and Hubble types from the Hyperleda database except 
otherwise noted. S? denotes unclassified spiral galaxy. 

a From Table II in Marziani ct al. (1999), where UGC 3995A is referred to as UGC 3995B. 
b Full name SDSS J141824.74+250650.7. 
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TABLE 6 

HI Gas Properties of other Galaxies in Active Galaxy Field 
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Note. — Neighboring galaxies in Group I are those identified from HI maps to be tidally inter- 
acting with their corresponding active galaxies, in Group II are those identified as probably tidally 
interacting with their corresponding active galaxies, and in Group III as possible candidates for 
producing the observed HI disturbances in their corresponding active galaxies (see text). 



a Full name SDSS J141824.74+250650.7. 
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Fig. 1. — Upper panel: Optical image of Markarian 1 (Mrk 1) (active galaxy) and NGC 451 from the 2nd Digitized Sky Survey (DSS2). 
Lower panels: (Left) Contours of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity- 
weighted HI mean velocity (first moment). Mrk 1 lies at the center of the map, whereas NGC 451 is located south-east of center. In the 
zeroth moment map, contours are plotted at 5, 10, 15, 20, 25, 30, 35 x25 mjy beam -1 km s _1 (1.25 X 10 19 cm -2 ). In the first moment 
map, velocities are indicated by the scale wedge, and contours plotted at intervals of 25 km s _1 . The ellipse at the lower right corner of 
the lower panels indicates the half-power width of the synthesized beam, and has a size of 61" X 54". 
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Fig. 2.— HI channel maps of Mrk 1 and NGC 451. Contour levels are plotted at -6, -4, -3, 3, 4, 6, 8, 10, 13, 17 X 0.43 mjy beam" 1 
(1<t), which corresponds to a HI column density of 4.43 xlO 18 cm -2 . The central heliocentric velocity is shown for each channel. The cross 
marks the position of Mrk 1, and star the position of NGC 451. The synthesized beam is shown by the ellipse at the lower right corner of 
the top left panel. 
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Fig. 3. — Upper panel: Optical image of NGC 7679 and NGC 7682 (both active galaxies) from the DSS2. Lower panels: (Left) Contours 
of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean velocity (first 
moment). NGC 7679 lies to the west and NGC 7682 to the east of center. In the zeroth moment map, contours are plotted at 1, 5, 10, 
20, 30, 40, 50, 70, 90 x20 mjy beam" 1 km s" 1 (9.12 X 10 18 cm~ 2 ). In the first moment map, velocities are indicated by the scale wedge, 
and contours plotted at intervals of 25 km s~ x . The ellipse at the lower left corner of the lower panels indicates the half-power width of 
the synthesized beam, and has a size of 68" X 53". 
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Fig. 4.— HI channel maps of NGC 7679 and NGC 7682. Contour levels are plotted at 2, 3, 5, 7, 9, 11, 15, 20, 25 X 0.45 mjy beam" 1 
(1<t), which corresponds to a HI column density of 4.23 xlO 18 cm -2 . The central heliocentric velocity is shown for each channel. The cross 
marks the position of NGC 7679, and star the position of NGC 7682. The synthesized beam is shown by the ellipse ar the lower left corner 
of the top left panel. 
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Fig. 5. — Upper panel: Optical image of ARK 539 (active galaxy) and PGC 2370360 from the DSS2. Lower panels: (Left) Contours 
of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean velocity (first 
moment). ARK 539 lies at the center of the map, whereas PGC 2370360 is located south-east of center. In the zeroth moment map, 
contours are plotted at 1, 5, 10, 15, 20, 30, 40, 50 X 15 mjy beam -1 km s _1 (7.86 X 10 18 cm -2 ). In the first moment map, velocities 
are indicated by the scale wedge, and contours plotted at intervals of 25 km s _1 . The ellipse at the lower left corner of the lower panels 
indicates the half-power width of the synthesized beam, and has a size of 58" X 54". 
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Fig. 6.— HI channel maps of ARK 539 and PGC 2370360. Contour levels are plotted at 2, 3, 4, 6, 10, 15, 20, 28 X 0.38 mjy beam" 1 
(1<t), which corresponds to a HI column density of 4.10 xlO 18 cm~ 2 . The central heliocentric velocity is shown for each channel. The 
cross marks the position of ARK 539, and star the position of PGC 2370360. The synthesized beam is shown by the ellipse at the lower 
left corner of the top left panel. 
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Fig. 7.— Upper panel: Optical image of NGC 5548 (active galaxy), SDSS J141824.74+250650.7, and PGC 1725892 from the DSS2. Lower 
panels: (Left) Contours of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI 
mean velocity (first moment). NGC 5548 lies at the center of the map, whereas SDSS J141824. 74+250650. 7 is located south-east and PGC 
1725892 south-west of center. In the zeroth moment map, contours are plotted at 1, 5, 10, 15, 20, 30, 40, 50 X 35 mjy beam -1 km s -1 
(1.94 X 10 19 cm -2 ). In the first moment map, velocities are indicated by the scale wedge, and contours plotted at intervals of 25 km s -1 . 
The ellipse at the lower left corner of the lower panels indicates the half-power width of the synthesized beam, and has a size of 56" X 53". 
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Fig. 8.— HI channel maps of NGC 5548 (active galaxy), SDSS J141824.74+250650.7, and PGC 1725892. Contour levels are plotted at 
2, 3, 5, 7, 9, 11 xO.45 mjy beam -1 (lc), which corresponds to a HI column density of 5.18 X 10 18 cm -2 . The central heliocentric velocity 
is shown for each channel. The cross marks the position of NGC 5548, star the position of SDSS J141824. 74+250650. 7, and triangle the 
position of PGC 1725892. The synthesized beam is shown by the ellipse at the lower left corner of the top left panel. 
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Fig. 9. — Upper panel: Optical image of Markarian 461 (Mrk 461) (active galaxy) and PGC 48870 from the DSS2. Lower panels: 
(Left) Contours of integrated HI intensity (zcroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean 
velocity (first moment). Mrk 461 lies at the center of the map, whereas PGC 48870 is located south-west of center. In the zeroth moment 
map, contours are plotted at 1, 5, 10, 20, 30, 40, 50 X 20 mjy beam -1 km s _1 (4.49 X 10 18 cm -2 ). In the first moment map, velocities 
are indicated by the scale wedge, and contours plotted at intervals of 25 km s _1 . The ellipse at the lower left corner of the lower panels 
indicates the half-power width of the synthesized beam, and has a size of 64" X 54". 
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Fig. 10.— HI channel maps of Mrk 461 and PGC 48870. Contour levels arc plotted at 2, 3, 5, 7, 9, 15 X 0.43 mjy beam -1 (la), which 
corresponds to a HI column density of 4.40 X 10 18 cm -2 . The central heliocentric velocity is shown for each channel. The cross marks the 
position of Mrk 461, and star the position of PGC 48870. The synthesized beam is shown by the ellipse at the lower left corner of the top 
left panel. 
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Fig. 11. — Upper panel: Optical image of NGC 841 (active galaxy) and NGC 834 from the DSS2. Lower panels: (Left) Contours 
of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean velocity (first 
moment). NGC 841 lies south-east of center, whereas NGC 834 is located north-west of center. In the zeroth moment map, contours are 
plotted at 1, 5, 10, 20, 30, 50, 70 X 30 mjy beam" 1 km s" 1 (1.39 X 10 19 cm" 2 ). In the first moment map, velocities are indicated by the 
scale wedge, and contours plotted at intervals of 25 km s _1 . The ellipse at the lower left corner of the lower panels indicates the half-power 
width of the synthesized beam, and has a size of 65" X 54". 
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Fig. 12.— HI channel maps of NGC 841. Contour levels are plotted at 2, 3, 4, 6, 10, 15, 20, 25, 30, 40 X 0.42 mjy beam" 1 (la), which 
corresponds to a HI column density of 4.01 X 10 18 cm -2 . The central heliocentric velocity is shown for each channel. The cross marks the 
position of NGC 841. The synthesized beam is shown by the ellipse at the lower left corner in the top left panel. 
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Fig. 13. — HI channel maps of NGC 834. Contour levels are plotted at 2, 3, 4, 6, 10 X 0.42 mjy beam 1 (lcr), which corresponds to 
a HI column density of 4.01 X 10 cm . The central heliocentric velocity is shown for each channel. The cross marks the position of 
NGC 834. The synthesized beam is shown by the ellipse at the lower left corner of the top left panel. 
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Fig. 14. — Upper panel: Optical image of Markarian 573 (Mrk 573) (active galaxy) and PGC 1226059 from the DSS2. Lower panels: 
(Left) Contours of integrated HI intensity (zcroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean 
velocity (first moment). Mrk 573 lies at the center of the map, whereas PGC 1226059 is located north-east of center. In the zeroth moment 
map, contours are plotted at 1, 5, 10, 15, 20, 25 X 20 mjy beam -1 km s _1 (9.62 X 10 18 cm -2 ). In the first moment map, velocities are 
indicated by the scale wedge, and contours plotted at intervals of 25 km s _1 . The ellipse at the lower left corner of the lower panels is the 
half-power width of the synthesized beam, and has a size of 67" X 51". 
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Fig. 15. — HI channel maps of Mrk 573 and PGC 1226059. Contour levels are plotted at 2, 3, 5, 7 X 0.5 mjy beam 1 (la), which 
corresponds to a HI column density of 4.95 X 10 18 cm -2 . The central heliocentric velocity is shown for each channel. The cross marks the 
position of Mrk 573, and star the position of PGC 1226059. The synthesized beam is shown by the ellipse at the lower left corner of the 
top left panel. 
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Fig. 16.— Upper panel: Optical image of NGC 7591 (active galaxy), PGC 214933, and HI 2318+065 from the DSS2. Lower panels: (Left) 
Contours of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean velocity 
(first moment). NGC 7591 lies at the center of the map, whereas PGC 214933 is located south-west and HI 2318+065 east of center. In 
the zeroth moment map, contours arc plotted at 1, 5, 10, 20 , 30, 40, 50, 70, 90, 110 x20 mjy beam -1 km s _1 (7.39 X 10 18 cm -2 ). In the 
first moment map, velocities arc indicated by the scale wedge, and contours plotted at intervals of 25 km s — 1 . The ellipse at the lower left 
corner of the lower panels indicates the half-power width of the synthesized beam, and has a size of 63" X 53". 
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Fig. 17.— HI channel maps of NGC 7591, PGC 214933, and HI 2318+065. Contour levels are plotted at 2, 3, 4, 6, 10, 20, 30, 40, 50, 
70 X 0.45 mjy beam -1 (lcr), which corresponds to a HI column density of 3.43 X 10 18 cm -2 . The central heliocentric velocity is shown for 
each channel. The bigger cross marks the position of NGC 7591, star the position of HI 2318+065, and the smaller cross the position of 
PGC 214993. The synthesized beam is shown by the ellipse at the lower left corner of the top left panel. 



Tidal Interactions among Scyfcrt Galaxies 



43 




Heliocentric velocity 
480Q 49QD BQQD 5100 



31 24 



* 22 

s 

CM 

I - 



z 
□ 



14 



20 kpc 



. 1 ■ 1 

■ * 

■ I 

• ■ 



IB 



16 - 




12 
02 28 45 



o 



_L 




o 



_L 



30 15 00 

RIGHT ASCENSION (J200C) 



27 45 



30 15 00 

RIGHT ASCENSION (J2000) 



27 45 



Fig. 18. — Upper panel: Optical image of Markarian 1040 (Mrk 1040) (active galaxy) and its companion galaxy PGC 212995 from the 
DSS2. Lower panels: (Left) Contours of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity- 
weighted HI mean velocity (first moment). Mrk 1040 is the larger galaxy at the center of the map, whereas PGC 212995 is the smaller 
galaxy located just to the north. In the zeroth moment map, contours are plotted at 1, 5, 10,20, 40, 80, 120, 160 X 30 mjy beam" 1 km s" 1 
(2.19 X 10 19 cm -2 ). In the first moment map, velocities are indicated by the scale wedge, and contours plotted at intervals of 25 km s — 1 . 
The ellipse at the lower left corner of the lower panels indicates the half-power width of the synthesized beam, and has a size of 50" X 54". 
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Fig. 19.— HI channel maps of Mrk 1040 and PGC 212995. Contour levels are plotted at 2, 3, 4, 6, 10, 20, 30, 50, 70 X 0.45 mjy beam" 1 
(1<t), which corresponds to a HI column density of 6.77 X 10 18 cm -2 . The central heliocentric velocity is shown for each channel. The cross 
marks the position of Mrk 1040. The synthesized beam is shown by the ellipse at the lower left corner of the top left panel. 
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Fig. 20. — Upper panel: Optical image of Markarian 341 (Mrk 341) (active galaxy) and NGC 169 from the DSS2. Lower panels: (Left) 
Contours of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean velocity 
(first moment). Mrk 341 is the smaller galaxy at the center of the map, and NGC 169 the larger galaxy located just to the north. In the 
zeroth moment map, contours are plotted at 1, 5, 10, 20, 30, 50, 70, 90, 110 X 20 mjy beam -1 km s _1 (1.24 X 10 19 cm -2 ). In the first 
moment map, velocities are indicated by the scale wedge, and contours plotted at intervals of 25 km s _1 . The ellipse at the lower left 
corner of the lower panels indicates the half-power width of the synthesized beam, and has a size of 54" X 49". 
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Fig. 21.— HI channel maps of Mrk 341 and NGC 169. Contour levels are plotted at 2, 3, 4, 6, 8, 10, 15, 20 X 0.5 mjy beam" 1 (la), 
which corresponds to a HI column density of 6.39 X 10 cm . The central heliocentric velocity is shown for each channel. The southern 
cross marks the position of Mrk 341, and the northern cross the position of NGC 169. The synthesized beam is shown by the ellipse at the 
lower left corner of the top left panel. 
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Fig. 22. — Upper panel: Optical image of UGC 1395 (active galaxy) and HI 0155+066 from the DSS2. Lower panels: (Left) Contours 
of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean velocity (first 
moment). UGC 1395 lies at the center of the map, whereas HI 0155+066 is located west of center. In the zeroth moment map, contours 
are plotted at 1, 5, 10, 20, 30, 40, 50 X 20 mjy beam -1 km s _1 (8.5 X 10 18 cm -2 ). In the first moment map, velocities are indicated by the 
scale wedge, and contours plotted at intervals of 25 km s — 1 . The ellipse at the lower left corner of the lower panels indicates the half-power 
width of the synthesized beam, and has a size of 73" X 53". 
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Fig. 23.— HI channel maps of UGC 1395 and HI 0155+066. Contour levels are plotted at 2, 3, 4, 6, 8, 10, 12 x 0.46 mjy beam" 1 (la), 
which corresponds to a HI column density of 4.03 X 10 18 cm~ 2 . The central heliocentric velocity is shown for each channel. The cross 
marks the position of UGC 1395, and star the position of HI 0155+066. The synthesized beam is shown by the ellipse at the lower left 
corner in the top left panel. 
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Fig. 24. — Upper panel: Optical image of Markarian 993 (Mrk 993) (active galaxy) and HI 0126+321 from the DSS2. Lower panels: 
(Left) Contours of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean 
velocity (first moment). Mrk 993 lies at the center of the map, whereas HI 0126+321 is located east of center. In the zeroth moment map, 
contours are plotted at 1, 5, 10, 20, 30, 50, 70, 90 X 20 mjy beam -1 km s~ 1 (3.92 X 10 18 cm -2 ). In the first moment map, velocities 
arc indicated by the scale wedge, and contours plotted at intervals of 25 km s . The ellipse at the lower left corner of the lower panels 
indicates the half-power width of the synthesized beam, and has a size of 79" X 53". 
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Fig. 25. — HI channel maps of Mrk 993. Contour levels are plotted at 2, 3, 4, 6, 10, 11, 20, 25 X 0.45 mjy beam 1 (lcr), which corresponds 
to a HI column density of 1.82 X 10 18 cm -2 . The central heliocentric velocity is shown for each channel. The cross marks the position of 
Mrk 993. The synthesized beam is shown by the ellipse at the lower left corner in the top left panel. 
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Fig. 26. — HI channel maps of HI 0126+321. Contour levels are plotted at 2, 4, 6, 10, 16 X 0.42 mjy beam 1 (lcr), which corresponds 
to a HI column density of 1.70 X 10 cm" . The central heliocentric velocity is shown for each channel. The cross marks the position of 
HI 0126+321. The synthesized beam is shown by the ellipse at the lower left corner in the top left panel. 
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Fig. 27. — Upper panel: Optical image of Markarian 3157 (Mrk 3157) (active galaxy) and HI 0446+185 from the DSS2. Lower panels: 
(Left) Contours of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean 
velocity (first moment). UGC 3157 lies at the center of the map, whereas HI 0446+185 is located north-west of center. In the zeroth 
moment map, contours arc plotted at 1, 5, 10, 15, 20, 30, 40 X 30 mjy beam -1 km s _1 (1.7 X 10 19 cm -2 ). In the first moment map, 
velocities are indicated by the scale wedge, and contours plotted at intervals of 25 km s _1 . The ellipse at the lower left corner of the lower 
panels indicates the half-power width of the synthesized beam, and has a size of 58" X 50". 
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Fig. 28.— HI channel maps of Mrk 3157 (active galaxy) and HI 0446+185. Contour levels are plotted at 2, 3, 4, 6, 10, 15, 20 X 
0.42 mjy beam -1 (lcr), which corresponds to a HI column density of 4.90 X 10 18 cm -2 . The central heliocentric velocity is shown for each 
channel. The cross marks the position of Mrk 3157, and star the position of HI 0446+185. The synthesized beam is shown by the ellipse 
at the lower left corner in the top left panel. 
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Fig. 29. — Upper panel: Optical image of Markarian 1419 (Mrk 1419) (active galaxy) from the DSS2. Lower panels: (Left) Contours 
of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean velocity (first 
moment). In the zeroth moment map, contours are plotted at 1, 5, 10, 15, 20, 30, 40 X 20 mjy beam -1 km s _1 (8.14 X 10 18 cm -2 ). In 
the first moment map, velocities are indicated by the scale wedge, and contours plotted at intervals of 25 km s -1 . The ellipse at the lower 
left corner of the lower panels indicates the half-power width of the synthesized beam, and has a size of 66" X 61". 
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Fig. 30. — HI channel maps of Mrk 1419 (active galaxy). Contour levels are plotted at 2, 3, 5, 7, 9 X 0.43 mjy beam 1 (lcr), which 
corresponds to a HI column density of 3.61 X 10 18 cm -2 . The central heliocentric velocity is shown for each channel. The cross marks the 
position of Mrk 1419. The synthesized beam is shown by the ellipse at the lower left corner in the top left panel. 
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Fig. 31. — Upper panel: Optical image of IRAS 14082+1347 (active galaxy) from the DSS2. Lower panels: (Left) Contours of integrated 
HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity- weighted HI mean velocity (first moment). In the 
zeroth moment map, contours are plotted at 1, 3, 5, 7 X 10 mjy beam -1 km s _1 (5.08 X 10 18 cm -2 ). In the first moment map, velocities 
are indicated by the scale wedge, and contours plotted at intervals of 25 km s _1 . The ellipse at the lower left corner of the lower panels 
indicates the half-power width of the synthesized beam, and has a size of 61" X 53". 
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Fig. 32. — HI channel maps of IRAS 14082+1347. Contour levels are plotted at 2, 3, 4 X 0.43 mjy beam 1 (la), which corresponds to 
a HI column density of 4.50 X 10 cm . The central heliocentric velocity is shown for each channel. The cross marks the position of 
IRAS 14082+1347. The synthesized beam is shown by the ellipse at the lower left corner in the top left panel. 
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Fig. 33. — Upper panel: Optical image of UGC 3995 (active galaxy) and UGC 3995A from the DSS2. Lower panels: (Left) Contours 
of integrated HI intensity (zcroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean velocity (first 
moment). UGC 3995 lies to the east, whereas UGC 3995A to the west of the center. In the zeroth moment map, contours are plotted at 1, 
5, 10, 20, 30, 50, 70 X 40 mjy beam -1 km s -1 (1.22 X 10 19 cm -2 ). In the first moment map, velocities are indicated by the scale wedge, 
and contours plotted at intervals of 25 km s _1 . The ellipse at the lower left corner of the lower panels indicates the half-power width of 
the synthesized beam, and has a size of 110" X 49". 
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Fig. 34.— HI channel maps of UGC 3995 and UGC 3995A. Contour levels are plotted at 3, 5, 7, 9, 11, 15, 20, 25, 30 xO.62 mjy beam" 1 
(1<t), which corresponds to a HI column density of 3.90 xlO 18 cm~ 2 . The central heliocentric velocity is shown for each channel. The 
cross marks the position of UGC 3995, and star the position of UGC 3995A. The synthesized beam is shown by the ellipse at the lower left 
corner in the top left panel. 
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Fig. 35. — Upper panel: Optical image of NGC 1167 (active galaxy) from the DSS2. Lower panels: (Left) Contours of integrated HI 
intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean velocity (first moment). NGC 1167 
lies at the center of the map, whereas HI 0302+352 is located to the north-east. In the zeroth moment map, contours are plotted at 1, 
5, 10, 15, 20, 30, 40, 50, 60 X 30 mjy beam -1 km s _1 (1.45 X 10 19 cm -2 ). In the first moment map, velocities are indicated by the scale 
wedge, and contours plotted at intervals of 25 km s _1 . The ellipse at the lower left corner of the lower panels indicates the half-power 
width of the synthesized beam, and has a size of 65" X 52". 
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Fig. 36.— HI channel maps of NGC 1167 and HI 0302+352. Contour levels are plotted at -9, -7, -5, -3, 3, 5, 7, 9, 11 15, 20, 
30 X 0.6 mjy beam - 1 (lcr), which corresponds to a HI column density of 5.97 X 10 18 cm" 2 . The central heliocentric velocity is shown for 
each channel. The cross marks the position of NGC 1167, and star the position of HI 0302+352. The synthesized beam is shown by the 
ellipse at the lower left corner in the top left panel. 
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Fig. 37. — Moment maps of IC 1749, detected in the same field as the target object UGC 1395. (Left) Contours of integrated HI 
intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity- weighted HI mean velocity (first moment). In the 
zeroth moment map, contours are plotted at 1, 5, 10, 15, 20, 25 X 35 mjy beam -1 km s -1 (1.49 X 10 19 cm -2 ). In the first moment map, 
velocities are indicated by the scale wedge, and contours plotted at intervals of 25 km s _1 . The ellipse at the lower left corner of each panel 
indicates the half-power width of the synthesized beam, and has a size of 73" X 53". 
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Fig. 38. — Moment maps of NGC 845, detected in the same field as the target object NGC 841. (Left) Contours of integrated HI 
intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity- weighted HI mean velocity (first moment). In the 
zeroth moment map, contours arc plotted at 1, 5, 10, 15 X 30 mjy beam -1 km s" 1 (1.39 X 10 19 cm -2 ). In the first moment map, velocities 



are indicated by the scale wedge, and contours plotted at intervals of 25 km s 
the half-power width of the synthesized beam, and has a size of 65" X 54". 
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Fig. 39. — Moment maps of PGC 3097128, detected in the same field as the target object UGC 3157. (Left) Contours of integrated HI 
intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity- weighted HI mean velocity (first moment). In the 
zeroth moment map, contours are plotted at 1, 3, 5, 10, 15, 20 X 35 mjy beam -1 km s -1 (1.83 X 10 19 cm -2 ). In the first moment map, 
velocities are indicated by the scale wedge, and contours plotted at intervals of 25 km s -1 . The ellipse at the lower left corner of each panel 
indicates the half-power width of the synthesized beam, and has a size of 58" X 54". 



HI 1346+341 



473 a 



Heliocentric velocity (km/s) 
4a 00 4020 4B4D 



3 



5 



3415 
14 
13 
12 
11 



□ 

(> 10 



hi 
□ 



09 



i 1 r 

20 kpc 



1 1 F 




13 40 4Q 





30 25 20 15 10 05 
RIGHT ASCENSION (J200O) 



00 



35 30 25 20 15 10 05 
RIGHT ASCENSION (J2000) 



00 



Fig. 40. — Moment maps of HI 1346+341, detected in the same field as the target object Markarian 461. (Left) Contours of integrated 
HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity- weighted HI mean velocity (first moment). In the 
zeroth moment map, contours are plotted at 1, 5, 10 X 8 mjy beam -1 km s -1 (3.80 X 10 18 cm -2 ). In the first moment map, velocities are 
indicated by the scale wedge, and contours plotted at intervals of 25 km s -1 . The ellipse at the lower left corner of each panel indicates 
the half-power width of the synthesized beam, and has a size of 64" X 54". 
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Fig. 41.— Moment maps of PGC 2369294 and HI 1828+503, detected in the same field as the target object ARK 539. (Left) Contours 
of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean velocity (first 
moment). PGC 2369294 lies at the center of the map, whereas HI 1828+503 is located north-west of center. In the zeroth moment map, 



contours are plotted at 1, 3, 5, 10, 15, 20 x35 mjy beam 1 km s 
by the scale wedge, and contours plotted at intervals of 25 km s _1 . 
width of the synthesized beam, and has a size of 58" X 55". 
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Fig. 42. — Moment maps of HI 2318+064, detected in the same field as the target object NGC 7591. (Left) Contours of integrated HI 
intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity- weighted HI mean velocity (first moment). In the 
zeroth moment map, contours are plotted at 1, 5, 10 X 20 mjy beam -1 km s _1 (7.39 X 10 18 cm -2 ). In the first moment map, velocities 
are indicated by the scale wedge, and contours plotted at intervals of 25 km s _1 . The ellipse at the lower left corner of each panel indicates 
the half-power width of the synthesized beam, and has a size of 63" X 53". 
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Tidal Interactions among Seyfert Galaxies 
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Fig. 43.— Moment maps of UGC 12628, detected in the same field as the target objects NGC 7679 and NGC 7682 . (Left) Contours 
of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean velocity (first 



moment). In the zeroth moment map, contours arc plotted at 1, 5, 10, 20, 30, 40 x40 mjy beam 1 km s 1 (1.82 X 10 19 



2 ). In the 



first moment map, velocities are indicated by the scale wedge, and contours plotted at intervals of 25 km s 1 . The ellipse at the lower left 
corner of each panel indicates the half-power width of the synthesized beam, and has a size of 68" X 53". 
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Fig. 44. — (Upper) Plot of number (solid line) and cumulative fraction (dashed line) of Seyfert galaxies in our ensemble sample with 
(candidate) interacting neighboring galaxies as a function of their projected separation up to 200 kpc. The data has been divided into bins 
of width 10 kpc. A vertical dotted line is drawn at 190 kpc, within which the galaxies are all mapped with a nearly uniform sensitivity 
(at worst a difference of ~20%). (Lower) Plot of the difference in radial velocity between the Sefyert and all detected neighboring galaxies 
within 210 kpc. Crosses denote those in Group I, triangles those in Group II, asterisks those in Group III, and boxes the remainder not 
identified to be interacting with their corresponding Seyfert galaxy. Horizontal dashed lines are drawn at ±100 km s _1 to help guide the 
reader. 



